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Abstract 



Analytic calculations of the Lamb shift represent a considerable challenge due to the size and the 
complexity of the expressions that occur in intermediate steps. In the current work, we present a method 
for the treatment of the bound-state self-energy in higher orders. Advantage is taken of computer algebra 
systems. The method is applied to the 2P-states, and the (one-loop) self-energy is calculated up to the 
order of a(Za) 6 . The calculation leads to improved predictions for the in atomic hydrogen. 

In particular, the focus of the current work is the calculation of the so-called Ago-coefficient with improved 
accuracy. Ago-coefficients are calculated with an accuracy of less than 1Hz in frequency units, 

A 60 (2P 1/2 ) = -0.99891(1) and A 60 (2P 3/2 ) = -0.50337(1). 

An account of the contributions to the Lamb shift of 2P states is given [status of 1996 when this thesis 
was completed]. Bound-state quantum electrodynamics combines the complexity of modern quantum 
field theory, including renormalization, augmented by the structure of bound states and the inherent 
relativistic corrections with the experimental possibilities of modern high-resolution laser spectroscopy. 



This thesis was submitted to the University of Munich on 10 February 1996, and final examinations 
took place on 8 May 1996. We attempt to provide a more complete account of the calculations whose 
results were previously presented in [U. D. Jentschura and K. Pachucki, e-print physics/0011008, Phys. 
Rev. A 54, 1853 (1996)]. The current version is not identical to the master thesis submitted to the 
University of Munich; it contains hypcrrefcrences and updates. Part of the results are of relevance for 
a recent investigation of radiative corrections induced by local potentials [U. D. Jentschura, e-print 
hep-ph/0305066, J. Phys. A 36, L229 (2003)]. 

The thesis is in German; this English abstract is provided for completeness. Final results for the Lamb 
shift of 2P states as reported in Ch. 6 represent the status of corrections known by 1996. 



This thesis is the first in the traditional threefold sequence "diploma(mastcr)-dissertation(PhD) - 
habilitation" that is being followed in Continental Europe. Copies of the author's PhD thesis are available 
as U. D. Jentschura, "Quantum Electrodynamic Radiative Corrections in Bound Systems (Drcsdncr 
Forschungen: Theoretische Physik, Band 2)" , w.e.b. Univcrsitatsverlag, Dresden, 1999 (ISBN: 3-933592- 
65-8), 225 pages (address of w.e.b. Publishers: BergstraBe 78, 01069 Dresden, Germany). Copies of 
the habilitation thesis (title: "Quantum Electrodynamic Bound-State Calculations and Large-Order 
Perturbation Theory") as well as copies of the master thesis in its 1996 version are available via e-mail or 
"snail mail" (jentschura@physik.uni-freiburg.de or U. D. Jentschura, Physikalisches Institut, Universitat 
Freiburg, Hermann-Herder-Strafie 3, 79104 Freiburg im Breisgau, Germany). 
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Zusammenfassung 



Analytischc Rcchnungcn zur Lamb-Vcrschicbung wcrdcn durch Groflc und Komplcxitat dcr auftrctcndcn 
Terme erheblich crschwcrt. Es wird in dieser Arbeit eine Methodc vorgcstcllt, mit wclchcr dcr 
Selbstenergiebeitrag zur Lamb-Vcrschicbung wasserstoffartiger Systemc in hohercr Ordnung analytisch 
berechnet werden kann. Durch den Einsatz eines symbolischen Computcrprogramms gclingt cs, die 
Rechenschritte zu vcreinheitlichcn und zu beschleunigen. 

Die entwickelte Methode wird auf die 2P-Zustande angewandt und der (Ein-Schleifen-) 
Selbstenergiebeitrag bis zur Ordnung a(Za) e berechnet. Die Rechnung zeichnet sich gegeniiber 
fruhcrcn numerischen Rechnungen durch deutlich geringere Fehlergrenzen aus. 

Insbcsondere ist es das Ziel dieser Arbeit, fur die 2P-Zustande den sogenannten Ago-Kocffizienten mit 
verbesserter Genauigkcit zu berechnen. Fur die 2P-Zustande ist dieser KocfBzicnt bisher nur durch 
Extrapolation einer numerischen Rechnung bekannt [P. J. Mohr. Phys. Rev. A 26, 2338 (1982)], 

A 60 (2P 1/2 ) = -0.8(3) und A 60 (2P 3/2 ) = -0.5(3). 

In dieser Arbeit werden die Ago-Kocfhzientcn mit einer Genauigkcit von < 1 Hz berechnet, 

A 60 (2P 1/2 ) = -0.99891(1) und A 60 (2P 3/2 ) = -0.50337(1). 

Dadurch kann die Lamb-Verschiebung fur die 2P-Zustande mit verbesserter Genauigkeit angegeben 
werden (bisher lag die theoretische Unsicherheit bei etwa 3 kHz). Es ergeben sich folgende neue 
theoretische Werte fur die Lamb-Verschiebung der 2P-Zustande, 

<5P Lamb (2P 1/2 ) = -12836.0(3) kHz 

sowie 

5P Lamb (2P 3/2 ) = 12517.5(3) kHz, 

jeweils mit einer theoretischen Unsicherheit von 300 Hz. Fur die Fcinstrukturaufspaltung der 2P-Zustandc 
wird als theoretischer Wert 

SE fs = P(2P 3/2 ) - P(2P 1/2 ) = 10969043.2(6) kHz 

erhalten. Ein Vergleich der Voraussagen der Quantenelektrodynamik mit dem Experiment 
(hochstaufloscndc Spcktroskopie) kann zur Ubcrpriifung dcr Quantcnfcldthcoricn diencn. 
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Kapitel 1 

Einleitung 



"Quantum Electrodynamics (QED) occupies a unique position in the hierarchy of theoretical physics. As 
the basic theory of the electromagnetic interaction, it provides the foundation for the quantum mechanics 
of atoms and molecules as well as condensed matter physics, not to mention their numerous applications 
in chemistry, biology, etc. At the same time, it is the first theory that has been applied successfully to 
high energy phenomena. " 
T. Kinoshita in [5] 



1.1 Die QED und die Lamb-Verschiebung 

Die Quantcnclcktrodynamik beschrcibt die Wcchsclwirkung des elektromagnctischen Feldes (des Lichtcs) 
mit geladenen Teilchen. Ausgangspunkt fur die Entwicklung der QED waren unter anderem experimentell 
nachgewiesene geringe Abweichungen von den Voraussagen der Diracschen Theorie des Wasserstoff atoms. 

Es sei daher an dieser Stelle an die Diracsche Formcl fiir die Encrgieniveaus des Wasserstoffatoms erinnert. 
Nach der Diracschen Theorie liegen die Niveaus der gebundenen Zustande gerade bei 

E n ,j = m e f(n,j) (1.1) 

mit 

/(n,i)=(l + ^ Y 1 ". (1.2) 

V (n - H + 1/2) + + V2) 2 " {Zaffl 

Eine gute Moglichkcit, die Diracschen Encrgieniveaus genauer zu untersuchen, liegt in der Entwicklung 
der Formel 1.2 in eine Taylorreihe in der Fcinstrukturkonstantcn a. In crster Naherung gilt 

E n ,j =m e - (Za) 2 m e /(2n 2 ) + 0(a 4 ), (1.3) 

d.h. bis zur Ordnung a 2 ist E n j von j unabhangig. Wir erkennen hier die Schrodingersche Formcl 
fiir die gebundenen Encrgieniveaus (E^ = —a 2 m e /(2n 2 )). Zu den Schrodingerschen Niveaus ist im 
relativistischen Sinne gerade die Ruhcnergic m e des Elcktrons addicrt worden. 

Zwischen den 2Pi/ 2 und 2P 3 / 2 -Nivcaus liegt nach der Diracschen Theorie cine Encrgicdiffercnz. In der 
Ordnung a ist diesc Differcnz das Fcinstrukturintervall 

4 

Cxm c 

E2P3/2 ~ E2P1/2 = Feinstrukturintervall = + 0(a ) fiir Z = 1. 

32 

In Einheiten der Frequenz ergibt sich daraus cine Differcnz von 

E 2P3 /2 - E 2P1/2 » 10949 MHz. 
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Die in der Schrodingerschen Theorie vorhandene Entartung der 2P-Niveaus wird also in der Diracschen 
Theorie aufgehoben. 

Andere Entartungen bleiben jedoch erhalten. So sind die Energieniveaus der Diracschen Theorie bei 
vorgegebenem Gesamtdrehimpuls j (j = I + s) vom Bahndrehimpuls I und Spin s unabhangig. Dahcr 
sind zum Beispiel die 2S , 1 / 2 und 2P 1 / 2 Niveaus in der Dirac-Theorie energetisch entartet, denn sic tragcn 
beide die Hauptquantcnzahl n = 2 und die Gesamtdrchimpuls-Quantcnzahl j — 1/2. Es gilt also in der 
Dirac-Theorie 

E 2P1 / 2 - E 2S i /2 = (exakt) 

und somit 

E 2P3 / 2 - E 2S1 / 2 = E 2P3/2 - E 2P1/2 = Feinstrukturintervall. 

W. Lamb und R. Retherford stellten jedoch 1947 in einer spektroskopischen Messung fest, dafi zwischen 
dem 2Pi/ 2 und dem 2S , i/ 2 -Niveau im Widerspruch zur Diracschen Theorie ein Energieunterschied besteht 
[10]. Die Differenz betragt nach einer verbesserten Messung von Triebwasser, Dayhoff und Lamb [11] aus 
dem Jahr 1953 

E 2Sl/2 - E 2Pl/2 = 1057.8(1) MHz. 
Dieser EfFckt wird hcute als "klassischc Lamb-Vcrschicbung" bezeichnct. 

Die Laserspcktroskopie crmoglicht heute spektroskopische Messungen an wasserstoffartigen Systcmcn 
mit einer relativen Gcnauigkeit von bis zu 10 -12 [12]. Der angestrebte Vcrglcich mit den Voraussagen der 
Quantcnclcktrodynamik crfordcrt Berechnungen zur Lamb-Verschicbung in hohcrcr Ordnung. 

Im Falle des Wasserstoffs kommen die Hauptbeitrage zur Abweichung von der Diracschen Theorie durch 
quantenelektrodynamische Effekte zustande. Einerseits sind kleine Abweichungen vom Coulomb-Gesetz 
bei sehr kleinen Abstanden dafiir verantwortlich (Vakuumpolarisation) , andererseits die Wechselwirkung 
des Elektrons mit dem virtuell vorhandenen eigenen Strahlungsfeld (Selbstenergie). Es gelang H. 
Bethe noch im Jahr der ersten Messung (1947), die Rechnungcn zu den Hauptbcitragcn der Lamb- 
Verschicbung (Vakuumpolarisation und Selbstenergie) in erster nichtverschwindender Ordnung (a(Za) 4 ) 
zu vervollstandigen [8] . 

1.2 Beitrage zur Lamb-Verschiebung 

Obwohl nicht alle Abweichungen von der Diracschen Theorie auf die Quantenelektrodynamik 
zuriickzufiihren sind, werden dennoch alle diese Abweichungen unter dem Begriff "Lamb-Verschiebung" 
zusammengefafit. Die Quantenelektrodynamik (QED), die Struktur des Kerns (Protons) und die voile 
relativistische Bchandlung des Wasserstoffs als Zwcikorpcrproblcm licfern aufier der Vakuumpolarisation 
und Selbstenergie noch verschiedenc andere Korrckturen zur Diracschen Theorie, die wir hier nach der 
Grofic ihrer Beitrage zum Wasscrstoffspcktrum auflistcn wollcn: 

• die Selbstenergie (self-energy), 

• die Vakuumpolarisation (vacuum polarization), 

• Ruckstofikorrekturen (recoil corrections), 

• Strahlungs-Riickstofikorrckturcn (radiative recoil corrections), 

• Zwci- und Drci-Schleifcn-Korrckturen zur Selbstenergie (two- and three-loop corrections), 

• Kerngrofienkorrektur (finite nuclear size effect). 

Es sci hicr der Klarheit halber erwahnt, dafi die Lamb-Verschicbung fur jeden Zustand des 
Wasscrstoffatoms cinzcln ausgewcrtet wird. Durch die Differenz der Verschiebungen der einzelenen 
Niveaus ergibt sich dann der Energieunterschied, der experimentell nachgewiesen wird. 

Man spricht bei alien hier erwahnten Effekten von (quantenelektrodynamischcn) Korrekturen zur 
Diracschen Theorie, weil die Effekte klcin sind und daher storungsthcoretisch behandelt werden konncn. 
Dies erleichtert die Rechnung bzw. macht die Auswertung der auftretenden Integrate iiberhaupt erst 
moglich. 
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1.3 Gegenstand der vorliegenden Arbeit 



''Historically, the one-photon radiative corrections (with all orders of the Coulomb interaction) have been 
considered in most detail, and the work required to attain the present accuracy has extended over four 
decades. " 
J. Sapirstein in [G] 



Das Thema der vorliegenden Arbeit ist der Ein-Schleifcn-Selbstenergiebeitrag zur Lamb-Verschiebung. 
Dieser Effekt kommt durch den Austausch eines virtuellen Photons zustande. Er stcllt den Hauptbeitrag 
zur Lamb-Verschiebung in Wasserstoff dar. Es soil eine Methode (efw-Methode) vorgestellt werden, mit 
der dieser Bcitrag in hoherer Ordnung analytisch berechnet werden kann (siehe Kapitel 3). Die Methode 
fiihrt zu einer gewissen Vercinfachung dcs sogenannten Nicdrigencrgic-Anteils. Der Hauptbeitrag wird von 
den relativistischen Korrekturen getrennt, und die Rcchnung kann in Tcilbcitrage aufgcspalten werden. 
Mit dieser Methode konnen fchlende Koeffizicntcn hoherer Ordnung berechnet werden. 

Die Methode wird in dieser Arbeit zur Berechnung der (Ein-Schleifen-) Selbstenergie der 2P-Zustande 
in der Ordnung a(Za) e angewandt. Unter Vernachlassigung der (geringen) Abhangigkeit von der 
reduzierten Masse des Systems, die spater miihelos wieder eingesetzt werden kann (siehe [6]), ist der 
Selbstenergiebeitrag zur Lamb-Verschiebung gegeben durch 

_ a m ( Za ) 4 jp 

O-t^SE TYl - — t, 

7r nr 

wobei sich fur den dimcnsionslosen P-Faktor die folgende Entwicklung in Potenzen von a ergibt 

F = F(Z)=A iQ + A 41 ln((Za)- 2 ) + A 5Q (Za)+ (1.4) 



(Za) 2 [A 6Q + A 61 \n((Za)- 2 ) + A 62 \n 2 ((Za)- 2 ) + O(Za)] . 

Die erste Ziffcr in den Indizcs der Koefhzientcn gibt die Potenz von (Za) an, die zweite Ziffer die Potenz 
des Logarithmus. Fiir P-Zustande verschwinden einige der Koeffizienten, insbesondere gilt A^i = A50 = 
A62 = [6]. Einige andcre Koeffizienten sind fiir die 2P-Zustande bekannt, 

A 4 o(2P 1/2 ) = -i-^lnfc (2P), 

bH 1/2) 18Q 

A4o(2P3 /2 ) = ^-^lnfco(2P), 
2Q 

A 61 (2P 3/2 ) = Vq 

lnfc (2P) = -0.0300167089(1). 

Fiir die ^go-Koeffizienten der 2P-Zustande existierten bis jctzt keine genauen Werte. Es bestcht jedoch 
cine Moglichkcit, die Werte dieser Koeffizienten mit Hilfc von numerischen Rcchnungcn abzuschatzen. 
Numerischc Rcchnungcn cxistieren fiir F(Z) lediglich fiir groBc Kcrnladungszahlen Z ', d.h. fiir schwerc 
wasscrstoffartige Ionen. Fiir kleine Z konvergiert die numerische Rechnung sehr langsam. P. Mohr [14] 
konnte in seiner Arbeit trotz grofieren Aufwands Ergebnisse nur fiir Z > 5 angeben. Durch eine 
Ausgleichsrechnung mit einer geeigneten Funktion und Extrapolation in den Bcrcich fiir kleine Z la8t 
sich dennoch eine Abschatzung fiir die j4go-Koeffizicntcn angeben. Die zur Ausgleichsrechnung benutztc 
Funktion lautet: 

G SE {Z) = A 60 + (Za)A 7Q 
Unter Venachlassigung des fiir Z = 1 kleinen Beitrags 

{Za)A 70 



sowie 



mit dem Bethe-Logarithmus [6] 
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ist dann 

G SE (Z = 1) w A 6Q . 

Die auf diese Art und Weise erhaltenen numerischcn Rcsultate fiir Aqq lauten: 

A 60 (2P 1/2 ) « G SE (2P 1/2 ) = -0.8(3) (1.5) 

sowie 

^6o(2P 3 / 2 ) « G SE (2P 1/2 ) = -0.5(3). (1.6) 

In dieser Arbeit sollen alle Koeffizienten, die in der Ordnung a{Zaf zur Lamb-Verschiebung 
beitragen, analytisch berechnet werden. Hierbei ist das Wort "analytisch" so zu verstehen, dafi die 
Rechnung im Hochenergie-Anteil (siehe Abschnitt 3.7) vollkommen analytisch, im Niedrigenergie-Anteil 
(siehe Abschnitt 3.8) nur bis auf verblcibcnde cindimcnsionalc Intcgrationen analytisch crfolgt. Dicsc 
(cindimcnsionalcn) Intcgrationen konvcrgicrcn sehr schncll und konncn mit groficr Gcnauigkcit ausgcfuhrt 
werden, da die Intcgrandcn glattc Funktionen ohnc Singularitatcn sind (Vcrfahrcn: Gaufischc Quadratur) . 

Die Ergcbnisse fiihren auf ncuc Wcrtc fiir die Lamb-Vcrschicbung der 2P-Zustande und einen neuen 
Wert fiir die Fcinstrukturaufspaltung. Wegen der Lange der Terme werden die meisten Zwischcnschrittc 
in den Rechnungen nicht explizit aufgelistet. Die GroBe der Ausdriicke iibersteigt bei den Integranden 
der Matrixelemente im Hochenergie-Anteil und bei den Zwischenschritten im Niedrigenergie-Anteil 
oft die Grenze von 1000 Termen. Selbst manche (Teil-) Ergebnisse sind so lang, daB nur deren 
Struktur (funktionalc Abhangigkeiten) angegeben wird. Es werden statt einer expliziten Angabe der 
Zwischenergebnisse die entscheidenden konzeptuellen Schritte des Verfahrens dargestellt. 

1.4 Motivation fiir die vorliegende Arbeit 

Die Quantenelektrodynamik ist bereits experimentell sehr gut bestatigt. Experimentelle Uberpriifungen 
der Quantenelektrodynamik haben besonderc Bcdcutung, denn die Thcoric ist Modcllthcorie fiir 
das Standardmodcll. Da cs nicht viclc Wahlmoglichkcitcn bei der Konstruktion renormierbarcr 
Quantcnfcldthcoricn gibt, wiirdc cine Abwcichung von den Vorhersagcn der QED zur Entwicklung ncucr 
Thcoricn fiihren musscn. 

Eine weitere Motivation fiir die aufwendigen Rechnungen crgibt sich aus der Prazision hcutigcr 
spektroskopischer Messungen. Die Genauigkeit der theoretischen Lamb-Verschiebung der iS-Zustande wird 
heute durch die Ladungsverteilung im Atomkern (den quadratischen Kernradius) beschrankt. Es besteht 
also die Hoffnung, durch neue Rechnungen und Experimente einige Erkenntnisse iiber Eigcnschaftcn der 
Elementartcilchcn zu gewinnen. 

Der in dieser Arbeit erhaltene genaue Wert fiir die Lamb-Vcrschicbung der 2P-Zustandc kann zum Beispiel 
als Rcfcrcnzpunkt zur Ubcrprufung der Mcssung der 2«S f -Lamb-Vcrschicbung diencn. 

1.5 Computer- Algebra Systeme 

Der Einsatz von Computer Algebra Programmen in der Theoretischen Physik hat es ermoglicht, 
umfangrcichc Bcrcchnungcn mit vielen tausend Termen in den Zwischenschritten in viel kiirzerer Zeit 
und sehr viel zuverlassigcr als bisher durchzufiihrcn. Symbolischc Umformungen werden vom Rcchncr 
vorgenommcn. Quantcnfcldthcorien bieten sich als Anwendungsgcbiet fiir computergestiitzte symbolischc 
Rechnungen in gewisser Weise an, denn storungstheoretische Rechnungen werden in hohercn Ordnungcn 
recht umfangreich. 

Auch in dieser Arbeit wird fiir die Berechnungen ein Computer Algebra System verwendet, das Programm 
Mathematica der Firma Wolfram Research. Dieses System zeichnet sich im Vergleich zu anderen 
Computer Algebra Systemen durch eine sehr flexible Programmiersprache aus, die sowohl regelbasierte 
wie auch prozedurale und funktionalc Elemente enthalt. Die Programmiersprache Mathematica ist auch 
objektorientiert: Mit den mathematischen Objcktcn konnen Regeln assoziiert werden. Auf diese Weise 
lasscn sich zum Beispiel Intcgrationsrcgcln definicrcn, die auch Ausdriicke mit mehreren tausend Termen 
im Intcgrandcn zuvcrlassig und schncll intcgrieren. Mit dicscm Programm konncn Matrixelemente der 
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Diracschen Wellenfunktion des 2P\/2 oder des 2P 3 / 2 Zustands automatisiert berechnet werden. Dabci 
entstehen Intcgranden mit mehreren tausend Termen. 

Es wurde ein weiteres Algebra-Programm verfafit, welches Kommutatoren von Operatoren automatisch 
auswertet. Mit Hilfe dieses Programms wird die Foldy-Wouthuysen Transformation des Dirac- 
Hamiltonians automatisch ausgewcrtet (mit bckanntcm Ergcbnis, siehe [3]). Das Algebra-Programm wird 
verwendet, um die Foldy-Wouthuysen Transformation des Operators a l e lk ' r = 7 7*e r auszuwerten 
(siehe Abschnitt 3.4). Das Ergcbnis diescr Transformation ist neu und ein zusatzlichcs Rcsultat dieser 
Arbeit. 
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Kapitel 2 

Theorie der Selbstenergie 



2.1 Konventionen 

Es folgt hicr cine Listc der am haufigsten verwendcten Symbole und Konventionen. 

• Es werden naturliche Einheiten verwandt, in denen h = c — 1 ist und somit nur eine einzige 
physikalische Dimension verbleibt. Als verbleibende Dimension wahlen wir in dieser Arbeit die 
Lange. Die Dimension der Lange wollen wir mit A e bezeichnen. 

• Griechische Indizes durchlaufen die Werte von ... 3, lateinische Indizes die Werte 1 ... 3. 

• Indizes von kontravarianten Vektoren werden als Supcrskript plaziert. Zum Bcispiel ist k° — uj die 
0-te Komponente des Vierervektors k. 

• Werden griechische oder lateinische Indizes wicdcrholt, so gilt die Summcnkonvention, auch wenn 
beide Indizes oben oder bcide untcn stehen (Beispicl a ■ p — a l p l ). Die Summcnkonvention wird 
aufgehoben, wenn der Index mit einem Strich versehen ist (zum Beispiel bei n' anstatt n). 

• Die Spurbildung einer Matrix erfolgt vermittels der Tr-Operation. 



2.2 Heaviside-Lorentz-Einheiten 

Um von SI-Einheiten zu Heavisidc-Lorcntz-Einhciten zu gclangcn, gchen wir aus vom Coulomb-Gcsctz 

F = fco^- = — 

r z 47reo r z 

und wahlen die Einheiten so, dafi eo — 1 wird. Aufierdem wahlen wir die Einheiten der Lange und Zeit 
so, dafi gerade c = 1 gilt (das heifit insbesonderc, dafi Lange und Zeit die gleiche physikalische Dimension 
besitzen). Wir wahlen die Einheit der Energie so, dafi auch h = 1 ist. In Heaviside-Lorentz-Einheiten ist 
dann gerade 

2 e 2 A k e 2 
e = — — = 47T— — = Ana 
eohc he 

(weil €q = h = c = I). Die elektrische Ladung ist in diesen Einheiten eine dimensionslose Grofie. Das 
anziehende Coulomb-Potential zwischen zwei Elcmcntarladungen wird dann 

. 1 e 2 1 Aira a 

V(r) = = = . 

47reo r 47reo t r 

Die erste Maxwell- Gleichung wird zu 

div E = — p = p, 

so dafi sich die inhomogencn Maxwcll-Glcichungcn in diesen Einheiten als 

d^ v =f (2.1) 

schrciben. 
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2.3 Die Selbstenergiefunktion 

Die Selbstenergic dcs Elcktrons kommt durch die Wechsclwirkung des Elcktrons mit dem cigenen 
Strahlungsfcld zustandc. Der quantenclektrodynamischc Prozcfi der Sclbstwcchsclwirkung lafit sich durch 
das Aussenden und das Wiedereinfangen cines virtuellen Photons bcschrcibcn. Zunachst wollen wir uns 
mit der Selbstenergiefunktion des freien Elektrons auseinandcrsetzen, die durch das Feynman-Diagramm 
2.1 beschrieben ist. Im Falle eines gebundenen Elektrons miissen wir den "einfach durchgezogenen " 




Abbildung 2.1: Sclbstcnergie des freien Elektrons 

freien Feynman-Propagator durch den vollcn Propagator des Elektrons, welches sich unter dem Einflufi 
des Zentralfeldes fortbewegt, ersetzen. Dies wird durch den Graphen 2.2 veranschaulicht, wobei die 
doppelt durchgezogene Linie den vollen Feynman-Propagator fur das Elektron (unter dem Einflufi dcs 
Coulombfeldes) darstellt: Zunachst wird die Selbstenergie des freien Elektrons besprochen. Gemafi den 




Abbildung 2.2: Selbstenergic dcs gebundenen Elcktrons 
Feynman-Regeln der QED ordnen wir diesem Graphen die folgende Amplitude zu. Wir schreiben 

• — ie7 M bzw. — iej u fur die zwei Vertizcs, 

• iDftv^k) fur den Photonen-Propagator zwischen den Vertizes, 

• iS F {p) fur den Feynman-Propagator dcs Elcktrons zwischen den Vertizes, 

5 und wir integrieren iiber den freien unbeobachteten Parameter, den Vierer-Impuls des Photons k. Wir 
konnen die Modifikation des Feynman-Propagators durch die Ein-Schleifen-Korrektur, die wir mit £(p) 
bezeichnen wollen, sofort hinschreiben [1]: 

f d A k 

-iE(p) = J j—^ iD^k) (-ie 7 M ) iS F (p) (-iej u ) 

oder umgcformt: 

m = ie?J ^D, v (k)rS F ( P )^. (2.2) 

Das Problem der "kleinen" Modifikation S zum Feynman-Propagator liegt nun darin, dafi S(p) unendlich 
grofi ist. Eine Abzahlung der fc-Potenzen ergibt 

d A k = dkO(k 3 ), 
D^(k) = Oik- 2 ), 
S F = Oik- 1 ), 
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dafi das Integral 2.2 linear divergiert. Dieses Problem wird durch Renormierung und Regularisierung 
gelost. Dazu wird der durch die Selbstwechselwirkung modifizierte Propagator des Elektrons 

Sp = Sf + Sp E(p) Sf + Sf E(p) Sf E(p) Sf + ■ ■ ■ 

gcnaucr untersucht. Wir konnen diese Rcihe formal aufsummieren zu: 

S F (P) = 1 y( r- . • (2.3) 

p — to — Zj(p) + it 

Beachtet man, dafi Sp(p) auf Spinor-Wcllcnfvmktionen wirkt, die zumindest naherungsweise auf der 
Massenschale liegen, so erscheint die folgendc Reihcnentwicklung um den Punkt p 1 = m sinnvoll: 

E(p) = A + B{i> - m) + E ron (p)(/f - to) 2 . 

Wir beachten, dafi E(p) und somit auch A, B und £ ren (p) von der Ordnung a sind. Wir approximicrcn 
mit dieser Vorgabc E(p) folgcndcrmaficn (Tcrme in hoherer Ordnung in a wcrdcn wcggclasscn): 

Slp ® = i> - m - A- B{j> - to) - (J - TO) 2 E rcn (p) ^ 2 ' 4 ) 



{p* - to - A)(l - B)(l -{p 1 - TO)E ren (p)) 
1 + B 

df-m-A)(l-(jJ-m)Eren(p)) 
(l+5)(l + (^-TO)E ren (p)) 

^ — TO — A 

^ (1 + B) 
p" — m — A 

wobei im letzten Schritt ausgenutzt wurde, dafi wir uns in der Nahe der Massenschale p 1 = m befinden. 
Nun erfolgt die Renormierung. Die Konstante 1+B wird in den Faktoren —ie-f^ der Elektron-Vertizes, die 
an den Selbstenergiegraphen angrenzen, wie folgt absorbiert. Man definiert die Renormicrungskonstantc 
Z2 als 

Z 2 = 1 + B (2.5) 
und die physikalisch observable Ladung des Elektrons als 

e'rt = Z 2 \fzle. 

Die ^-Renormierung ist also eine weitere Ladungsrcnormierung zusatzlich zur Z3-Ladungsrcnormierung 
( Vakuumpolarisation) . 

Das Elektron verhalt sich unter dem Einflufi seines eigenen Strahlungsfeldes offenbar so, als hatte es 
die Gesamtmasse to + A. Da man aber den Einflufi dieses Strahlungsfeldes im Experiment nicht einfach 
abschalten kann, so mufi man A als Massenrenormierung ansehen. Es ist also 

A = 5m. 

Rcgularisiert man den Photonenpropagator gemafi der Vorschrift von Pauli und Villars, 



1 



1 



1 



k 2 ~M 2 ' 



so ergibt sich 



3to 

47T 





1" 








+ 2 



(2.6) 



(2.7) 
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Dieser Ausdruck ist einerseits von der Ordnung a, also "klein" , andererseits ist er divergent in M, jedoch 
nur logarithmisch. Man absorbiert ihn in der renormierten Elcktroncnmasse. 

Um die Renormicrungskonstante B zu berechnen, miissen wir den Photonenpropagator weiter 
modifizieren gemafi 

-I -> 1 - 1 (2 81 

fc 2 fc 2 - fi 2 k 2 - M 2 ' y ' ' 

wobei der urspriinglichc Propagator im Limes fi — > 0, M — > oo crhalten wird. Wir erhalten dann 

Z 2 = 1 + B = 1 - 

Das Ergcbnis fiir Z 2 ist sowohl fiir fi — > als audi fiir M — > oo logarithmisch divergent. AuBcrdem ist 
der Wert fiir Z 2 nicht eichinvariant und hangt von der fiktiven Elcktronmasse fi ab. Wir werden aber im 
nachstcn Abschnitt sehen, dafi sich die nicht cichinvariante Z2-Renormicrung der Selbstencrgic gerade mit 
der ebenfalls nicht eichinvarianten Zi-Rcnormicrung der Vertexkorrektur weghcbt und daher die Thcoric 
konsistent bleibt. 

Die renormicrte Sclbstenergicfunktion crgibt sich dann als 

£ ren (p) = S(p) -5m- (Z 2 - m). (2.10) 



5 ln (^) +ln (^) + 5 



2.4 Die Vert exfunkt ion 

Durch die Selbstwcchsclwirkung cincs Teilchens mit seincm cigencn Strahlungsfcld vcrandcrt sich die 
Amplitude, die einem Wcchsclwirkungs- Vertex zugeordnet wird (Vertexkorrektur). Das entsprcchcnde 
Feynman-Diagramm ist gcgeben durch Abb. 2.3. Wir miissen daher 



Abbildung 2.3: Vertexkorrektur 

-iej^ -> -167^ - ieTf+ip,^) (2-11) 

ersetzen. Dabei sind p und p' die Viererimpulse dcs Elcktrons vor bzw. nach der Wechsclwirkung. Fiir 
die Vert exfunkt ion kann mit den Fcynman-Rcgeln leicht der folgcndc Ausdruck abgeleitet werden [1]: 

r„(P, P ') = -KfJ ^l^jz^'^j^-)- P-u) 

Auch die Vertexfunktion Y ^ mufi renormicrt werden, denn eine Abzahlung der fc-Potenzcn im Intcgrandcn 
crgibt cine logarithmische Divergenz des Integrals fiir T^. Aus der Entwicklung dcs Propagators 

(2.13) 



p 1 '— ft — vtl p 1 — ft — m + — p") 

1 1 



-(/- i)^— u (2.14) 

— ft — m p— ft — m p— ft — m 
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ergibt sich unter Abzahlung der fc-Potenzen der (divergente) Hauptbeitrag zur Vertexfunktion gerade fur 

p' = P- 

Daher kann T M geschrieben werden als 

r^{p',p) = T^( P ,p) + r^{p',p) (2.15) 

mit der renormierten Vertexfunktion T^(p',p). Man kann nachweisen, dafi der divergente Anteil der 
Vertexfunktion gerade proportional zu 7 M ist: 

T fl (p,p)=^L (2.16) 

Daher mufi eine weitere Ladungsrenormierung gemafi 

-ie7 M -> -ie-f^ - icy^L - ieT*(p,p') = -ie(l + L)T^{p,p) (2.17) 

vorgenommen werden, so dafi sich die drcifach rcnormiertc Elcktronladung zu 

4 = Z^e' R (2.18) 

ergibt mit 

Zy l = l+L odcr Z 1 = l-L (2.19) 

(in niedrigster Ordnung in a) und 

e' R = Z 2V %e (2.20) 

(^-Renormierung aufgrund der Selbstenergie, Z3-Renormierung aufgrund der Vakuumpolarisation). Mit 
der Ward-Identitat 

dpv 



wird unter Ausnutzung der Entwicklung von in Potenzen von {p 1 — m) 

£(p) = A + B(i> - m) + 0({p' - m) 2 ) 

die wichtigc Bcziehung 

= r M (p,p) = = -Bl» + 0{i> - m) (2.22) 

und somit 

L = -B (2.23) 

hergeleitet. Daraus folgt, dafi sich der bei der ^-Renormierung zur Elektronladung addierte Anteil und 
der bei der Vertexkorrektur-Renormierung addierte Term gerade herausheben: 

e' R = (1 + L + B)yfzle = (1 + (L - L))e = V^e. (2.24) 

Die nicht eichinvariantcn Z\- und ^-Renormierungen der Elektronladung heben sich gcgenseitig auf. 

2.5 Die Selbstenergie in einem gebundenen Zustand 

Der aufgrund der Selbstenergie modifizierte Feynman-Propagator ist 

Sp = Sf + Sf^tctiSf + SirE ren iS'i?S r eii<S'F + ■ ■ • (2.25) 

(das Supcrskript / bei 5*^ soil andeuten, dafi hicr die Intcraktion mit dem eigenen Strahlungsfeld 
beriicksichtigt wurde). 

Der Feynman-Propagator fur die Dirac-Glcichung lautct im Ortsraum, nach der Zeit fouriertransformicrt: 

S F {r u r 2 ,E) =3 — — — — . 2.26 

n hi — hi n (l — It) 
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Dabei erstreckt sich die verallgemeinerte Summe sowohl iiber den diskreten Teil des Spektrums 
(normalc Summation) wie auch iiber den kontinuierlichen Teil des Spektrums (Integration). Fiir diese 
verallgemeinerten Summen hat J. Schwinger das S-Symbol eingefuhrt. 

Die Pole des Propagators liegen gerade bei den Energie-Eigenwerten E n der Dirac Gleichung. 
Beriicksichtigt man die leichte Modifikation des Propagators durch S ren , so verschieben sich die Pole 
geringfiigig. Die neuen Positionen E n = E n + <j n (E n ) entsprechen dann den durch die Selbstencrgie 
vcrschobenen Energicniveaus. Wir beobachten (iiber n wird nicht summiert!) 

$»|7°WllU 



d 3 n / d 3 r 2 S ll+iflWrnfoW+faWnfa) 



(2.27) 



mn E-E„ 



E — E n 



und 



(^|7 ^ F S rcn 5 F 7V n ) = (^|7 5 F7 7 Src„^F7 |^> (2.28) 

= SS (^n\7 S Fl °\^)(^ Icn \ij m )(tP m \S Fl \ij n ) 



m I 



c C 1 - 1 

^9 p — w Snl 1 Eren I ^™)~? — 

ml hi — ±j„ 



Daher ist in erster Ordnung 



(E-E n y 



E — E n {E-E n y E-E n - (l/j n | ErenlV'n) 

und man kann in erster Ordnung die Energieverschiebung aufgrund der Selbstenergie als 

^se = (^„|£ rcn |V>„> (2.30) 

ansetzen. 

Aus dieser Gleichung fiir die Energieverschiebung eines gebundenen Zustands und der Gleichung 2.2 fiir 
die Selbstenergiefunktion ergibt sich der Selbstenergiebeitrag zur Lamb-Verschiebung im Zustand \tf>) als 

SEbe = (#e 2 J -0- A D^{k)rS v F { P ) Y -Ren|^) (2.31) 



wobei 

S % = 2 U 1 oT7 (2-32) 

den durch das Coulombpotcntial zusatzlich modifizierten Feynman-Propagator darstellt. Dieser 
Propagator wird auch als Dirac-Coulomb-Propagator bczcichnct. 

Ren steht in 2.31 fiir all diejenigen Ausdriicke, die bei der Renormicrung abgezogen werden miissen. 
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2.6 Renormierung der Selbstenergie 



Bei der Renormierung der Sclbstcnergicfunktion S(p) traten sowohl die Elcktron- Renormierung als 
auch die Massenrenormierung 5m auf. Die Glcichungcn 2.7 und 2.9 fur die Elcktron-Rcnormierung 
Z2 und die Massenrenormierung 5m warcn crhalten wordcn, indcm zunachst das Integral fur S(p) 
vor der Integration rcgularisicrt wurdc. Das Ergcbnis sind dann die logarithmisch divcrgcntcn 
Massenrenormierugs und Elektronrenormierungskonstanten. Es wurde dabei folgende Reihenfolgc 
eingehalten: Regularisierung, Integration, Renormierung. 

In praktischen Rcchnungen ist es jedoch haufig giinstiger, zunachst einen divergenten Ausdruck fur einen 
ProzeB zu berechnen und dann vom Integranden Terme abzuziehen, die den Renormierungen entsprechen. 
SchlicBlich rcgularisicrt man den Integranden und intcgricrt. Rcihcnfolge hierbci: Renormierung, 
Regularisierung, Integration. Diesc Reihenfolgc wollcn wir auch in dicscr Arbeit cinhalten und iibcrlegcn, 
welchen Ausdruck wir in dicsem Fall vom Integranden abzichen miissen. 

Eine Renormierung und Regularisierung der Selbstenergie des Elcktrons in cincm gebundenen Zustand 

5E SE = $|E(p)M (2.33) 

ist nur fur hohe Energien des ausgetauschten virtucllcn Photons (im ultraviolctten Limes) notwendig. 
Die Selbstenergicfunktion des gebundenen Elcktrons wird durch das elektrostatische Zentralkraftfclds 
modifizicrt, denn der frcic Propagator fur das Elcktron muB durch 

1 1 



j}— ji — m fl— ft — m — 7°y 

ersetzt werden. Fiir den Hochencrgic-Antcil der Selbstenergie ( "ultraviolctter Antcil") wahlen wir die 
Feynman-Eichung des Photoncnpropagators. Wir schreiben diesen Anteil, den wir mit Eh bezeichnen 
wollen, folgendcrmafien auf: 

Wir legen uns hierbei nicht fest, iiber welchen Bereich sich die Integration erstreckt. Dieser Punkt wird 
im Abschnitt iiber den Hochcncrgic- Anteil in Kapitel 3.7 naher besprochen. In der Formcl fiir Eh steht 
rcn fiir dicjenigen Ausdriickc, die bei der Renormierung der Selbstenergie abgezogen werden miissen. 

Wir definiercn das Matrixclcment 

P = till 1 * . u \ ^77mI^>- 

p— f — m — 7 U V 

AuBerdem fiihren wir fiir den Nenner im Feynman-Propagator die Abkiirzung 

D=jS-)i-m (2.35) 

ein. 

Fiir hohe Energien ui des virtuellen Photons ist es moglich und sinnvoll, den Dirac-Coulomb-Propagator 
in Potenzen des Potentials V zu entwickcln. Dabei wird ausgenutzt, dafi fiir gebundene Zustande im 
Coulombpotential V — — ~ von der Ordnung a 2 ist (r als Abstand vom Kern ist von der Ordnung 
a-Bohr = I /{am)). 

Die Entwicklung des Dirac-Coulomb-Propagators geschieht mit Hilfe der Formel 



X -Y X XX XXX 
Wir setzen X =f— ft — m = D und Y = j°V und erhalten als Resultat der Entwicklung 4 Terme, den 
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0- Vertex, 1- Vertex, 2- Vertex und 3- Vertex- Anteil. 

1 



(2.36) 



yl-k'-^V-m D-j°V 
1 

= D + 

l/vl + 

D D 

i/yl/yl + 
D D D 

I/yl 7 <yl 7 °V- + 
D D D D 

0(V 4 ) 

Hohcre Terme tragen nicht bei, da V 3 bereits von der Ordnung a 6 ist und wir nur bis zu dicscr Ordnung 
rcchncn wollcn. 

Diese Entwicklung entspricht einer Aufglicdcrung der Selbstcnergic dcs gcbundcnen Elcktrons in die 
Summe aus 4 einzelnen Feynman-Diagrammen gemafi der folgenden Abb. 2.4 Der erste dieser Graphcn 



Abbildung 2.4: Entwicklung der Selbstcnergic in Termen des Coulombfcldcs 

(der 0- Vertex- Anteil) entspricht gerade der Selbstcnergic des freien Elcktrons, der zweite dieser Graphcn 
(der 1-Vertex- Anteil) gerade der Vertexkorrektur zur einfachen Interaktion des Elcktrons am Coulombfcld. 
Die weiteren Graphen aus Abbildung 2.4 tragen keine Divergenzen, wic cine Abzahlung der fc-Potcnzen 
crgibt. 

Dahcr miissen wir bei der Rcnormicrung der Gcsamt-Selbstenergie sowohl die Rcnormicrungen aufgrund 
dcs 0- Vertex- Anteils wie auch aufgrund des 1- Vertex- Antcils mitcinbeziehen. 

Zunachst wird der 0- Vertex- Anteil behandclt, in wclchcm wir mit dem Propagator des freien Elcktrons 
rcchncn. Fur unscrc Rcchnungcn benotigen wir cine Integraldarstcllung des Masscn-Gcgcntcrms. Als 
Ansatz wahlcn wir 



Sm = —ie 



(2tt)4 



1 1 



k 2 k 2 _ M 2 



2 ^ m + ^ (2 37) 



Die Giiltigkeit dieser Darstellung sollte explizit nachgewiesen werden. Man kann den Ausdruck 2.37 
explizit berechnen und erhalt als Ergebnis Gl. 2.7. Im folgenden soli em Uberblick iiber die wichtigsten 
Rechentricks gegeben werden, die zur direkten Berechnung des Ausdrucks 2.37 angewandt werden. In der 
Rechnung wird zweckmafiig zu Hilfsvariablcn ct\ und cti iibergegangen, die gemafi 



k 2 + ie 

und 



/>oo 

/ dai exp(zo;i(fc 2 + ie)) 
Jo 



-i I dai cxp(iai(fc 2 — 2mk ■ t)) 



n 



k' 2 — 2muj 

eingefiihrt werden mit dem zeitartigen Einhcitsvcktor t 

t= (1,0,0,0) 
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Der Faktor uj kann als 



. d 



cxp(iz ■ k) 



z=0 



geschrieben werden. Mit der verallgemeinerten GauBschcn Integralformel 

d 4 k 



exp(i(afc + b ■ k)) = 



(47r) 2 a 2 



exp(-i67(4a)) 





1" 


( m 2 ) 






+ 2 



kann dann die 4-dimensionale ^-Integration ausgefuhrt werden. Dabci ist natiirlich der 
Photoncnpropagator gemafi Pauli-Villars zu regularisicren. Das Ergcbnis der Integration ist 

. 3m 
om = a 

47T 

Der Abzug des Massen-Gegenterms nach der Integration ist also aquivalent der Subtraktion des Ausdrucks 
2.37 vor der Integration. 

Es ist eine Besonderheit fur gebundene Zustande, da8 der mit der Elektronrenormierung Z 2 verbundene 
Term nicht beitragt und auch nicht abgezogen zu werden braucht. Unter Vernachlassigung des Massen- 
Gegenterms ist der divergente Anteil des O-Vertex-Diagramms 



mit 



B = Z 2 -1. 



(2.38) 



Der divergente Bcitrag aufgrund des 1- Vertex- Anteils ist, da zwischen den Feynman-Propagatorcn gcrade 
der Ausdruck j°V steht, gegeben durch 



mit 



und somit 

mit L = Z^ 1 - 1. Mit der Ward-Identitat 



Di = $\T°(p,p)V\1>) 
Di = $\fLV\il>) 



(2.39) 

(2.40) 
(2.41) 



d 



(2.42) 



lafit sich nun (siehe Abschnitt 2.4) die Beziehung 2.23 

L=-B 

zeigen. Daher gilt: 

D + D 1 = (i,\B(p'-m)\iP) + $\ 1 LV\i>) 

= (^|(Z 2 -l)(^-m-7°V)|V) =0, 

weil die Dirac-Gleichung mit Feld 

(y-m-7°V)|V') = 

erfiillt. Es ist daher lediglich notig, den Massen-Gegenterm Sm abzuziehen, und es ergibt sich in Feynman- 
Eichung: 

d A k 1 ,-, „ 1 



E H = -ie" 



(2tt) 4 k 2 



# - m - j°V 



(2.43) 
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Kapitel 3 

Die epsilon-fw-Methode 



3.1 Allgemeines 

Die Rechenmethode fur die Selbstenergie ist eine modifizierte Version der in [19] angewandtcn 
Rechenmethode (e-Methode). Die Methode stiitzt sich auf eine Trennung dcr Rcchnung in einen Hoch- und 
eincn Nicdrigcncrgicanteil (in Bczug auf die Encrgic dcs virtucllcn Photons). Diesc Auftcilung war bcreits 
Gegenstand dcr urspriinglichcn Rechnungen zur Selbstenergie gebundener Zustandc. H. Bcthc bcrcchncte 
in [8] den Niedrigenergie-Anteil fur die Selbstenergie in der Ordnung a(Za) 4 . Die hier vorgestellte 
Rechnung geht zwei Ordnungen in Za holier. Das bedeutet, daB in dcr Wcllcnfunktion und im Propagator 
relativistische Korrekturen zu der Betheschen Rechnung mitgenommen werden miissen. 

Die e-Methode wird in dieser Arbeit modifiziert. Im Niedrigenergie-Anteil wird eine Foldy-Wouthuysen- 
(FW)-Transformation cingcfuhrt. Dadurch crfolgt eine klare Trennung des Hauptbeitrages von den 
rclativistischen Korrekturen. AuBcrdcm kann mit dem nichtrclativistischcn Propagator fur das gebundene 
Elektron gcrcchnet werden. Der nichtrclativistischc Propagator ist in geschlossener Form sowohl im 
Ortsraum wic auch im Impulsraum bekannt. Die Rechnung lafit sich aufgrund dcr FW- Transformation in 
Einzelschritte aufspaltcn. die nachcinandcr abgcarbeitct werden konnen. Dadurch wird eine zusatzliche 
Vereinfachung erzielt (die Ausdriicke in den Zwischenschritten iibersteigen nur sehr selten die Grenze von 
1000 Tcrmen). Im Hochenergie-Anteil wird die Rechnung fast vollstandig computergestutzt ausgefiihrt. 
Zunachst wird in Potenzen des Potentials V und des Impulsoperators p (und damit in (Za)) entwickelt. 
Die dabei entstehenden Terme enthalten Matrixelemente mit der Diracschen Wellenfunktion \ip). Die 
Matrixelemente werden mit symbolischen Programmen zur beschleunigten Berechnung von Integralen 
ausgewertet (wobei hier die Grenze von 1000 Termen im Integranden oft iiberschritten wird). Als Name 
fur diese computer-unterstutze Foldy-Wouthuysen-e-Mcthode wird "efw-Methode" vorgeschlagen. 



3.2 Einige Definitionen 



Zunachst sei an die 4 7-Matrizen 7^ (/i = 0, 1, 2, 3) mit ihren Kommutatorbeziehungen 



{7^ , 7^} = 2.9' 



(3.1) 



crinnert. Wir vcrwenden fur diese 4 x 4-Matrizen die Dirac-Darstcllung, 




(3.2) 



Fur die 7' 



5 -Matrix verwenden wir die Konvcntion 



7 5 = 7°7 1 7 2 7 3 - 



(3.3) 



Wir wollen eine Matrix (einen Operator) als gerade definicren, wenn cr von der Form 




(3.4) 
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ist (mit 2 x 2-Matrizen A und B), d. h. wenn er in einer Bispinor-Darstellung Diagonalgestalt hat. Wir 
wollen eine Matrix (den Operator) ungerade nennen, falls er von der Form 



I I 



ist. Es sind zum Beispiel die a 4 -Matrizen mit 



(It) ("» 

ungerade (die a-Matrizen sind nicht mit der Fcinstrukturkonstantcn a zu verwechscln). Auch d-p = a l p l 
ist ein ungerader Operator. Dcmgcgcmiber ist zum Beispiel die 7°-Matrix gcradc. 

Nun zur Nomcnklatur. Wir wollen den Schrodingcr-Hamiltonian fur das Wasserstoffatom definieren als 

Hs = ^ + V. (3.7) 

Er wirkt auf dem Raum der komplexwcrtigcn Wcllenfunktioncn. Der Dirac-Hamiltonian 

H D = a-p + fim + V (3.8) 

wirkt auf dem Raum der 4-komponentigcn Dirac-Spinoren. Nun sei zur Klarheit darauf hingewiesen, daB 
wir beziiglich der Entwicklung in a fur gebundene Zustande 



a = 


0(a) 


(3.9) 


V = 


0(a) 


(3.10) 


(3m = 


0(1) 


(3.11) 


V = 


0(a 2 ) 


(3.12) 


\k\ = 


0(a 2 ) 


(3.13) 



anschen konncn. Fur das Potential V zum Beispiel ergibt sich die Bcgrundung folgcndernmaficn: r als 
Abstand vom Kern ist von der Ordnung ciBohi = l/(am), und wir beziehen uns in dieser Arbeit auf die 
Massenskala des Elektrons to. Daher fiihren wir den Parameter p mit 

tmr (3-14) 



GtBohr 

ein und crhaltcn 

V = -± = 

r 

(siehe das Buch Bethe und Salpeter, [9]). 



1G 



3.3 Die gamma-Matrizen 



In dieser Arbeit werdcn haung die 16 T-Matrizen vcrwcndct, die cine Basis der 4 x 4-Matrizen bilden. Es 
sind dies (siehe zum Beispiel [2], Seite 54): 

r° = I (3.15) 
r 1 - 4 = 7" (3.16) 
r 5 - 7 = i 1 ° 1 l =ia l (3.17) 

r 8 - 10 = l -e l3k ~f 3 i k = S l (3.18) 

T 11 - 14 = 7 V (3.19) 

T 15 = 7 5 (3.20) 

Diese 16 Matrizen werden als r a -Matrizen bezeichnet, wobei der Index a von ... 15 lauft. Die T a 
Matrizen ergeben sich durch Anwenden der Metrik auf die 7-Matrizen, aus denen die T-Matrix 
konstruiert ist. Zum Beispiel ist Tq = T° und r$ = — T 5 . Auf den 4 x 4-Matrizen kann mit Hilfe der Spur 
ein Skalarprodukt konstruiert werden, 

^ s > = E(i Tr(r " A) ) (i Tr ( r « B ))- ( 3 - 21 ) 

Mit der Beziehung 

Tr(r Q r' 3 ) =A8 a > 3 (3.22) 

rechnet man leicht nach, dafi die T-Matrizen eine Orthonormalbasis beziiglich des oben definierten 
Skalarprodukts bilden: 

(r Q , TP) = S al3 . (3.23) 
Das bedeutet aber, dafi jede 4x4 Matrix A in den T-Matrizen entwickelt werden kann, 

A = ^ V? (V? , A) = £ (±Tr(r A)) . (3.24) 

Diese Beziehung ist zum Beispiel im Buch von Itzykson und Zuber (siehe [2]) zu linden. Die Entwicklung 
in die T-Matrizen kann auf ein Skalarproukt vermittels der Spurbildung zuriickgefiihrt werden. Fiir die 
Berechnung der Spuren Tr(r^Mj) wird ein symbolischcs Programm verwandt [26]. 



3.4 Die Foldy-Wouthuysen-Transformation 

Wir wollen an das allgemeine Prinzip erinnern, auf dem die Foldy-Wouthuysen-(FW)-Transformation 
bcruht . Gcgeben sei ein Hamiltonian fiir Spinoren , 

H = <D + £ + f3m, (3.25) 

mit O als ungeradem Operator und £ als geradem Operator (so kann zum Beispiel im Dirac-Hamiltonian 
O = a ■ p und £ = V identifiziert werden). H ist hermitesch und somit auch O und £: 

+ = und £+ = £. (3.26) 

Wir haben aufgrund der expliziten Gestalt der Matrix 7 die Eigcnschaft 

{(3 , O} = oder f30 = -0/3. (3.27) 
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Zusatzlich sei H nicht explizit zeitabhangig, d. h. 

dH 



t =0. (3.28) 



Da H nicht zeitabhangig ist, transformiert sich die "Schrodingergleichung' 



dt 

untcr unitaren Transformationcn U wie 



H\ijj} (3.29) 



i^p- = UHU + {U\^)) falls ^=0. (3.30) 

Dies konnen wir auch so auffassen, als wiirde der transformierte Hamiltonian H' = U HU + auf die 
transformierte Wellenfunktion ip' = Uip wirken. 

Es konnen folgende Aussagen gezeigt werden. 

1. Der durch 

U = e lS mit S = -0 — (3.31) 

2m 

gegebene Operator ist unitar. 

2. Im transformicrtcn Hamiltonian 

H' = UHU + (3.32) 
verschwinden ungerade Operatoren in niedrigster Ordnung in a. 



Fiir die Unitaritat von U reicht es zu zeigen, daB S hermitesch ist. Denn dann ist 

U+U = e~ iS+ e iS = l. (3.33) 
Die Hermitezitat von S rechnet man leicht nach: 

S + = i^ + =i^ = -i^=S. (3.34) 

Der transformierte Hamiltonian H' schreibt sich unter Bcnutzung der Campbell-Baker-Haussdorff- 
Idcntitat als 

l" 



H' = e lS He- lS = ^- n [5,i/] (3.35) 

ra=0 

mit n [S, H] als 7i-Kommutator von S und H (° [S, H] = H) . In erster Ordnung ist daher 

H' = H + i [S, H] + . . . , (3.36) 
wobei wir den Kommutator [S, H] wiederum in niedrigster Ordnung in a als 

[S, H] = [-0— ,0 + £ + 13m] (3.37) 
2m 

O 
Im 

= iO + ... (3.38) 

schreiben k5nnen. Somit 

H' = H + i[S, H]=0 + £ + f3m + i(iO) + ■■■ = £ + [3m + .. . (3.39) 

Wie wir sehen, ist im transformicrtcn Hamiltonian der ungerade Operator in niedrigster Ordnung O 
gcrade vcrschwunden. 
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Um audi ungerade Operatoren hoherer Ordnung in a verschwinden zu lassen, miissen zwei und mehr FW- 
Transformationcn hintcrcinandergeschaltet werden. Der dann die Transformation vcrmittclnde Operator 
UU'U" ... ist dann selbstverstandlich immcr noch unitar. Dieses Verfahrcn wollen wir jetzt auf den 
Dirac-Hamiltonian anwenden. Wir identifizieren 

<D = a-p sowie £ = V. (3.40) 

Es ist dann 

S = -0-^- = ~0%r-- (3.41) 
Ira Ira 

In unserem Bcispicl reichen zwei FW-Transformationcn aus, um auch in Ordnung a 4 ungerade Operatoren 
zum Verschwinden zu bringen. Fur die erste der FW-Transformationen erhalt man nach einer langeren 
Rechnung bis zur Ordnung a 4 (siehe etwa [3]): 

/ O 2 O 4 \ 1 O 3 



2 



2m 8ra 3 I 8m 2 2m 3m 



1 (a-v) 3 

<a-p 7 [a-p,V]} + ^[a-p,V]-^-. (3.42) 



8m 2 ' ' 2m ' 3m 2 

Fiir die zweite FW- Transformation ist daher 

S = - 2/3 2m = ~* 2m" { ' £] ~ 3W ) = "»2^ W a ■ P ' V] - I^T ) ■ ( 3 ' 43 ) 

und nach einer ebenfalls langeren Rechnung erhalt man fiir den Foldy-Wouthuysen-transformeierten 
Hamiltonian i?FW (siehe zum Beispiel [3]) 

Hfw = H' = 7 ° ( m + f- - <f£) +V+ ^S(f) + -?L-Z • L. (3.44) 
\ Ira 8 m J J Zm z Am^r^ 

Die Foldy-Wouthuysen-(FW)-Transformation ist also nichts weiter als die Hintereinanderausfiihrung 
zweier unitarer Transformationen U und U', die durch exp(iS) und exp(zS") vermittelt werden und 
die jeweils bis zur Ordnung a 4 mitgenommen werden. Es ist 

H FW = U' (U HdU^U'- 1 = {U'U)H D (U'U)- 1 . (3.45) 

An der Glcichung 3.44 ersieht man im iibrigen, wic sich die relativistischen Korrekturn zur 
Schrodingerschen Theoric (Spin-Bahn-Kopplung, Zitterbewegungsterm und Darwin- Term) aus der 
Diracschen Theorie ergeben. 

Was konnen wir iiber die transformicrte Wcllenfunktion U\il>) aussagcn? Die Funktion \ip) = U'U\ip) ist 
bis zur Ordnung a 4 Eigenzustand des FW-Hamiltonians. Der FW-Hamiltonian ist in dicscr Ordnung 
einc gcradc Matrix, die obere und unterc Komponenten nicht mischt. Der Hauptantcil zum Eigenwert ist 
durch den Bcitrag 

' m 
— m 



7°m 



gegeben. \tp) ist Eigenzustand zu £/fw mit einem Eigenwert f» m. Daher mufi \ip) bis zu einer dem 
FW-Hamiltonian entsprechenden Ordnung frei von unteren Komponenten sein. 

Diese Ordnung cntspricht in unserem Fall der fiihrenden Ordnung der Dirac- Wcllenfunktion und ihrer 
ersten relativistischen Korrektur. Wir konnen daher in den folgcnden Rcchnungcn stcts die unteren 
Komponenten des Spinors l^) vernachlassigen und 7 — > 1, £ — > a setzen. 

Auf \tp) angewandt, ist der FW-Hamiltonian also gleichzusetzen mit der oberen Bispinor-Komponcntc 
dieses Operators, die wir mit bezeichnen wollen. 
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Im folgenden wollen wir nur noch mit der Anwendung des FW-Hamiltonians auf die oberen Komponenten 
eines Spinors beschaftigen und lassen daher das Superskript up wieder weg. Wenn wir uns an die Definition 



-2 



des Schrodingcr-Hamiltonians Hs = + V crinncrn und als rclativistische Korrcktur dazu 

SH = ~ ITT + TriW + 1 — ITT^ ' 1 ( 3 - 47 ) 
8 m 3 2m 2 4m 2 r 3 

definieren, so schreibt sich der FW-Hamiltonian sehr kompakt als 

H FW = m + H s + SH. (3.48) 

Die Ruhmasse m kann aus dem obigen Hamiltonian wegtransformiert werden. In nulltcr Ordnung ist 
die nichtrelativistische Schrodinger-Wellenfunktion \<p) Eigenzustand des FW-Hamiltonians. Durch die 
relativistischen Korrekturen SH kommt eine Korrektur \5<f>) hinzu. Die Normierung 

(0 + 5(f>\4> + 5<t>) = 1 in 1. Ordnung (3.49) 
fuhrt (wenn man (<^|<S<^>) als reell annimmt) auf 

{<f>\5<i>) = 0. (3.50) 

Die Korrektur SE zum Energieeigenwert ergibt sich aus 

SE = (4>\5H\4>). (3.51) 

\4> + 8(f)) soil die transformicrtc Schrodingerglcichung (mit relativistischen Korrekturen) approximativ 
losen, d. h. 

(E s + SE) |0 + 50) = (H s + 5H) \<j> + 5<j>) in 1. Ordnung (3.52) 
Es ergibt sich die folgende Differentialgleichung fur \8<j>), 

(H - (4>\5H\<t>)) \4>) = (E S - H S ) \8<t>), (3.53) 

und nach weiterer Umformung 

1^) = iE8 l Hs y *m (3-54) 

mit 

G^Es - H S ) = {Es ^ Hs y ( 3 - 55 ) 

als reduzierter Greensfunktion. 

Das Ergcbnis der oben aufgefiihrten FW- Transformation des Dirac-Hamiltonians ist bekannt. Es soil nun 
die FW- Transformation der Ausdriickc 

y 3 = oP cxp(ifc • r) (3.56) 

berechnet werden (diese Ausdriicke sollen mit derselben unitaren Transformation transformiert werden 
wie der Hamiltonian). Wir erhalten 

VU = U> (UyW- 1 ) U'- 1 = (V U) y 3 (V U)-\ (3.57) 

wobei U und U' diejenigen unitaren Transformational sind, mit denen auch der FW-Hamiltonian 
transformiert wurde. Da k eine a 2 -Potenz tragt und r von der Ordnung a -1 ist, konnen wir y 3 in 
Potenzen von a entwickeln. 

Vl = a 3 , (3.58) 
y{ = a 3 (ik-r), (3.59) 

Vi = oj(-\(k-r)f. (3.60) 
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Wir nehmen jeweils die um a 2 hShere Korrektur zum yQ-Term noch mit. Nach einer langeren Rechmmg, 
die hier nicht vorgefiihrt werden soil, erhalt man als Ergcbnis fiir die FW- Transformation von y^: 



<™ = « , +7 u ^-^rj (3.61) 

Vi 

y { FW = a 3 (ik [« ■ P , * " A ~ 2ikJ + 2 (k ' »V) ( 3 - 63 ) 

und schlicfilich fur 



, i) 



^- (2 7 °((fc • r)V + [P J , (ft ' r) 2 } - 7V [a -p,(k- r) 2 )] 



Fiir die Berechnung der Lamb-Verschiebung konnen in den letzten drei Ergebnissen die ungeraden 
Operatoren weggelassen werden, d. h. die zu den a-Matrizen proportionalcn Tcrme werden vernachlassigt. 
Man erhalt nach weiteren Umformungcn als Bcitrage der geraden Operatoren 

d fw = 7° - ^f) - (?* (3.65) 
' \m zw 5 J 2m z r A V / 

sowie 

^W = -^7 (ftx^ + l/(ft.f)y ( 3. 66 ) 

und 

»U = (* • f ) 2 ^ + (* • f ) (* x g )' ■ (3 - 67) 

Da wir uns wie oben nur mit den oberen Komponenten der Spinoren beschaftigen, konnen wir wie beim 
FW-Hamiltonian 7 — > 1 und S — > er setzen (wir beschranken uns auf die linke obere 2 x 2-Submatrix). 
Wir erhalten dann 

^={i-^P 2 )-^ir^) 3 (3-68) 

sowie 

y[, FW = -tL (k x <tY + i- (ft • f) pi (3.69) 



2m V / m 
und 

1/1™ - ~ (k-r) 2 p> + ±(k-r)(kx of. (3.70) 
Unter Vernachlassigung ungcrader Operatoren und fiir die linke obere 2 x 2-Submatrix konnen wir dahcr 

y> FW = U' (Ua^e^U- 1 ) U^ 1 = d FW + y{ FW + d FW (3.71) 



setzen mit 



P 3 (, , ..ft A 1 (t ,\ 2 



- 77-^ -3 ( r x °y 



2 m 3 2 m 2 r 3 

— — (k • f \ (k X a\ — ( ft x 

2ml A J 2 m 
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3.5 Der P hot onen- Propagator in verschiedenen Eichungen 

Der Photonen-Propagator ist defimert als Losung der Glcichung 

A^y) = f d 4 xD^(y - x\f(x). (3.73) 



Im Impulsraum gilt daher 

A fl (k)=D^(k)f{k). (3.74) 
In Lorentz-Eichung lauten die Maxwell- Glcichungen gerade 

UA»{x)=f{x), (3.75) 

und im Impulsraum folgt daraus in Lorentz-Eichung 

UA^k) = -k 2 A^k) = j^x) = -k 2 D^(k)f(k), (3.76) 

so daB sich 

Dny = —^rrr in Lorcntz- oder Feynman-Eichung (3.77) 
k z 

ergibt. Wir konnen zum Propagator beliebige Ausdriicke der Form 

addicren, ohne dafi sich an den beobachtbarcn Groficn F^ u = d^A u — d y A^ etwas andert. Bei der Rechnung 
mufi man nur die Kontinuitat dcs Stromes 

= o 

ausnutzen und bedenken, daB im Impulsraum 

dfi * ~ik^ 

zu ersetzen ist. In der Coulomb-Eichung 

diA\x) = (3.78) 

lauten die Maxwell-Glcichungen 

UA^x) - d^dMx) = j^x) (3.79) 



oder 

Es ergibt sich nach einer kurzen Rechnung fiir die Coulomb-Eichung 



k 2 A^k) + k^k A (k) = j (k). (3.80) 



D 00 (k) - ^, (3.81) 
D 0l {k) = Ao = 0, 

Die Transformation zwischen Feynman- und Coulomb-Eichung wird durch die Funktionen 

fo = sowie fi = -^r (3.82) 
k k l 

vcrmittclt. 
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3.6 Die epsilon-Methode 



Der urspriinglichc Ausdruck fur den Sclbstcncrgicbcitrag sicht cinfach aus, 

5Ese = ie2 J wr^^ ^-^-ln-^v ^ - <*i*»w ■ ( 3 - 83 ) 

Bei einer vollkommen analytischen Rechnung entstehen jedoch sehr viele Terme. Dies erschwert die 
analytische Behandlung. Aufierdem ergeben sich sowohl infrarote als auch ultraviolette Divergenzen, 
so dafi eine Teilung des Integrationsbereichs sinnvoll erscheint. 

Ein weiterer Grund fur die Einfiihrung einer Trcnnung des Intcgrationsweges ist darin zu suchen, dafi 
im relativistischen Bereich {lo w to) der Dirac-Coulomb Propagator in Potenzen des Coulombfcldes V 
entwickelt werden kann. Im nichtrclativistischen Bereich jedoch (u) <C m) tragen alle V Potenzen des 
Propagators bei. 

Es wird daher in Bezug auf die Energie to des virtuellen Photons eine Teilung des Integrationsweges 
vorgenommen. Man teilt die Rechnung in zwei Bercichc ein, die durch den Parameter e <C m voneinander 
getrennt sind und wahlt fur die beiden Bereiche unterschiedliche Eichungen des Photonenpropagators. 
Das Zicl ist, fur beide Bereiche eine effiziente Berechnungsmcthode zu wahlen und dadurch die analytische 
Integration zu crleichtcrn. 

Im Ausdruck fur die Selbstenergie ist der Dirac-Coulomb Propagator 

nV _ 1 

F j- # - m - 7 °V 

cnthalten. Wegen 

{ip ist Eigcnzustand des Hamiltonians), konncn wir im Ausdruck fur die Selbstenergie den Dirac-Coulomb- 
Propagator durch 

$F — ¥ \ ^r-^ K77 = Sf(E^ - lu) 

7 U (A^ — uj) — 7 • p — m — 7 U V 
ersetzen. Der Propagator Sp(E^ — lo) hat, wie man an der Darstellung 2.26 



S%(E-u) =S 



n (E^ - lo) - E n (l - irf) 
crkcnnt, Pole bei 

{Eli, — lo) — E n {l — ir/) = odcr lo = E^ — E n (l — it]). 

Unter Vernachlassigung des diskreten Spektrums, welches fur unsere Rechnung kcinc Rolle spielt, da die 
erzeugten Pole oberhalb der reellen Achse liegen, gibt es Energieeigenwerte E n fur E n > m und E n < —m. 
Fiir unscrcn Fall ist aufierdem E^ < m, E^ f» m. Daher gibt es Pole als Funktion von lo gcrade bei 

8*(w) < 0,9f(w) = it] = E n > 

sowie 

a?(w) > 2m, SS(cj) = -ir] = E n < 0. 

Da also im ursprunglichcn Elcktronen-Propagator die Pole fiir die positiven Eigenwerte E n leicht unter 
die reelle Achse gedriickt werden und die Pole fiir negative Eigenwerte einen leicht positiven Imaginarteil 
besitzen, werden die Pole fiir lo > 2m unter die reelle Achse gedriickt, und die Pole fiir lo < bekommen 
einen leicht positiven Imaginarteil. 

Es ergeben sich auch Pole des Integranden aus dem Photonenpropagator. Der Photonen-Propagator ist 
in Feynman-Eichung 

9iw 9iiv 



ir] 
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Aus der Umformung 

1 _ 1 

lu 2 — k 2 + it] (lu — \k\ + irf)(u> + \k\ — irj) 
crhaltcn wir fiir die Positionen der Singularitaten: 

lu = \k\ — irj sowie w = — \k\ + if). 

Diese Pole liegen daher fiir 9? (a;) > leicht unterhalb der reellen Achse und fiir 5R(w) < leicht dariiber. 



Jm(w) 

Pole des diskreten Spektrums^ 

WXX X 


fr/ 

' < ^F / |C ,2m 


Re{io 







Abbildung 3.1: Intcgrationswcgc fiir die ^-Integration 

Die w-Integration erfolgt unter Beriicksichtigung der Lage der Polstellen entlang des Feynman-Weges Cf 
(Abb. 3.1). Das Matrixelement P im Integranden geht fiir groBe k wie fc _1 , der regularisierte Propagator 
wie k _i und d 4 k geht wie A: 3 (zusammen k~ 2 ). Daher konncn wir Cf durch cincn uncndlich groBen 
Halbkreis in der unteren komplexen Halbebene schliefien (Lange dieses Halbkreises wie fc 1 , daher geht 
das Halbkrcis-Intcgral fiir groBe k wie fc _1 ). 

Man verbiegt nun die Intcgrationskontour Cf zu der modifiziertcn Kontour Cl + Ch (dabei wcrdcn kcine 
Residuen ausgelassen) . Die Integrationskontour C = Cl + Ch wird wie in der Abbildung 3.1 angegeben 
in der vorliegenden Arbeit verwandt. Da C infinitesimal iiber bzw. unterhalb der reellen Achse verlauft, 
konnen wir sicherlich Cl durch Hinzunahme cincs Vcrbindungsweges bei 5R(aj) = e schliefien. 

Wir definiercn als Hochencrgic- und Niedrigcnergic-Antcile 

Jc H ( 2?r ) P~ fi-m- 7 U K 

und 

El = ie ' 2 I wh D ^ \ u 1 m?^- 

J Cl (27r) 4 p p- j/t — m- TV 

Rcgularisierung und Renormierung sind nur fiir den Hochcncrgic-Antcil notig. Die Summe der beiden 
Antcile ist der gesuchte Selbstcnergie-Beitrag zur Lamb-Shift, 

SEse = El + Eh ■ 

5Ese hangt natiirlich nicht von der Wahl des (kleinen) Trennungsparameters e ab, da durch diesen 
Parameter nur der Integrationsbcreich gctrennt wird. Wir haben es also mit zwei Funktioncn El und Eh 
zu tun, die in Abhangigkcit von der Fcinstrukturkonstanten a und dem Parameter e betrachtet wcrdcn, 
und deren Summe nicht von e abhangt: 

5E S e(o) = E L {a,e) + E H {a,e) 

Die Entwicklungskoeffizienten von SEse beziiglich a hangen nicht von e ab. Der Entwicklungskoeffizient 
von SE beziiglich a ist aber gerade die Summe der entsprechenden Entwicklungskoeffizienten von El und 
Eh- Wir cntwickeln daher El und Eh zuerst in a. Dann wird jeder Entwicklungskoeffizient von El und 
Eh beziiglich a nachtraglich in e entwickelt. Die Summe dieser Entwicklung(cn) beziiglich e darf aber nicht 
von e abhangen, es miissen sich alle e-Kocffizienten (bis auf den von e°) herausheben. Die Koeffizienten 
von e tragen daher zur Selbstcncrgie bei, hohere und niedrigere (divergentc) Koeffizienten heben sich 
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gerade heraus. Daher geniigt es, die Entwicklungskoeffizienten von El und Eh beziiglich a nachtraglich 
bis zur Ordnung e° zu entwickeln. Terme, die fur e — ► logarithmisch oder starker divergieren, heben sich 
gerade heraus. 

Dies ist auch ein wichtiger Test fur die Rechnung: Derm jeder Rechcnfchler bci El und Eh wirkt sich 
mit hohcr Wahrschcinlichkcit auf den Kocfnzicnten dcr ln(e) oder 1/e-Divcrgcnz aus. Das Herausheben 
der Divergenz (in der Summe von El und Eh) ist somit ein wichtiger Test dcr Rechnung. 

Es ist fur die Giiltigkcit dcr Methodc wichtig. dafi weder El noch Eh von dcr Eichung abhangen. Dazu 
rcicht es aus, die Eichinvarianz von El zu zeigen. Da SE$e = El + Eh nicht von der Eichung abhangt, 
crgibt sich daraus automatisch die Eichinvarianz von Eh- Es ist also zu zeigen, dafi durch Modifikationcn 
am Propagator der Form (siehe Abschnitt 3.5) 



kcinc Veranderung des e° Kocfnzicnten von El cintritt. Unter Benutzung dcr Bcziehung 

1 m = f^! +m+7 !!+ /; m ^°Uw 



J& — m — 7°y \ ji— j/i — m — "f°V ft — m — j°V t 

= (-1 + 0)^) = -!^) (3.84) 
konnen wir den eichabhangigen Anteil £;£ auso umformen: 

Er sc = -^ 2 J CL j^^u k »+M* ( 3 - 85 ) 

p— ft — m — 7 u v / 

Dieser Ausdruck verschwindet fur e — > und tragt somit nicht zur Lamb-Shift bci. Daher ist der 
entscheidende e° Koefhzient in El eichinvariant. 

Als Eichung wahlen wir fur den Hochcnergic- Anteil den gemafi Pauli-Villars regularisierten Propagator 
in Fcynman-Eichung (siehe Abschnitt 3.5) 



i f 



k 2 fc2_ M 2 

fur den Niedrigenergie- Anteil jedoch die Coulomb-Eichung (ebenfalls siehe 3.5). 



3.7 Der Hochenergie- Anteil 

Wir schrciben den Hochenergie- Anteil Eh mit der expliziten Form des regularisierten Photonen- 
Propagators (in Feynman-Eichung) und mit der Massenrenormierung 5m noch einmal auf, 



E H = -ie z 



1 



f,2 k 2 _ M 2 



1 



f- # - m - -y°V 



7/d 



'{ip\5m\ip). 



(3.86) 
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Um diesen Ausdruck bis zur Ordnung a(Za) 6 zu berechnen, ist es notig, das Matrixelement 



p = ^1^-5—5 ou-rM ( 3 - 87 ) 

p— fc — m — y'V 

bis Ordnung a 6 auszuwerten. Wir wollen fur die spateren Rechnungen (tp\ explizit als ("0I7 schreiben. 
Damit ist 

Wie bereits erwahnt, ist die Integrationskontour Ch nach unten durch e begrenzt (Infrarot-Rcgulator). 
Daher ist cs moglich und sinnvoll, den Dirac-Coulomb-Propagator in Potenzen des Potentials V zu 
entwickcln. Fiir gcbundene Zustande im Coulombpotential ist V = — |f von der Ordnung a 2 (siehe 
Abschnitt 3.2). Wir miissen daher das Matrixelement P bis zur Ordnung V 3 entwickeln. 

Um diese Entwicklung bcsondcrs ubcrsichtlich ausfiihrcn zu konncn, vercinbarcn wir fiir den Ncnncr im 
Feynman-Propagator die Abkurzung (siehe Abschnitt 2.4) 

D=$-#-m. (3.89) 

Die Entwicklung des Dirac-Coulomb-Propagators geschieht mit Hilfe der Formel 

1 =L + 1 7Y L + L Y L Y L + ... 



X -Y X XX XXX 

Wir setzen X =i>— ft — m = D und Y = 7°V und erhalten als Resultat der Entwicklung 4 Terme, den 
0- Vertex, 1- Vertex, 2- Vertex und 3- Vertex- Anteil. 



1 1 



Jt - ^oy _ m D _ ^oy 



(3.90) 



1 

= D + 

W- + 

D D 

Hohere Terme tragen nicht bei, da V 3 bereits von der Ordnung a e ist und wir nur bis zu dieser Ordnung 
rechnen wollen. Wir definieren dann als Matrizen 

Mo = 7 V (3-91) 

Mr = 7 V-lyvi 7M (3.92) 

M 2 = 7 vI 7 0yl/^I 7Ai (3.93) 

M 3 = / 7 ^I 7 °vi/V-i 7 vl7M- (3-94) 

Wir konnen somit das Matrixelement P entwickeln als 

P = P +P 1 +P 2 +P 3 (3.95) 



2G 



mit 

P = (V|7°7"^7mM = WW) (3-96) 

Pi = M7V^7°^7„M = WMiW (3-97) 

h = ^\l°l^l°V^ 1 V^ 1 ^) = ^\M 2 \^) (3.98) 

P 3 = (^^l^V^V^V^M = (^|M 3 |V). (3.99) 

Diese Matrixelemente (0- Vertex, 1-Vertex, 2- Vertex und 3- Vertex- Antcile) sind auszuwerten. Dicsc 
Auswertung erfolgt. ctwas vcrcinfacht beschricben, in den folgcnden Schrittcn: 



1. Entwicklung der M^-Matrizen in die 16 T-Matrizen. Gemafi der Gl. 3.24 gilt: 

Mi = ^a p T p (3.100) 



mit 

o /J = i'&(r /J M i ). 

Die Entwicklung in die T-Matrizen macht die Rechnung iibersichtlich und systematisiert in gewisser 
Weise das Vorgehen, was sich besonders beim Einsatz des Computers positiv bemerkbar macht. 

2. Entwicklung der a^-Kocfhzientcn in Potcnzcn der Feinstrukturkonstanten a. Die 
Entwicklungskocffizientcn ap sind rationale Ausdriickc in den Impulsen des Elektrons, im 
Viererimpuls des virtuellen Photons und im Potential V und kSnnen daher in Potenzen von a 
entwickelt werden (da jeder der Operatoren eine bestimmte Potenz von a tragt). Es gilt: 

6 

a p = ^2a J bj(a p ). (3.101) 
j=o 

Es wird also jeder a^-Kocffizient bis a 6 entwickelt. 

3. Mittelung uber die raumlichen Anteile des Photoncnimpulses. Uber die raumlichen Komponenten 
in den &j(a/3)-Koeffizienten von k kann gemittelt werden, da sich die Integration uber den gesamten 
R 3 erstreckt. 

bj(ap) = f(k, andere Arg.) — > f(\k\, andere Arg.) = bj(ap). (3.102) 
Dabci entsteht der Ausdruck b(ap), wclcher nur noch von |fc|, to und naturlich von p und V abhangt. 

4. Ersetzung der Matrixelemente. Es werden nun (fur jede T^-Matrix einzeln) die 6j(a/3)-Koeffizienten 
durch die entsprechenden Matrixelemente mit \ip) ersetzt, gemafi 

bj{a p ) -+ (iPlTP'bjMlii), (3.103) 
(kcinc Summation uber /?, daher /?', siehe Konventionen) . 

Es entstehen aus den vorangegangenen Schrittcn Bcitrage wie zum Bcispiel (fur die Matrix Mi, 
Entwicklungskocffizient von 7 ): 

&4(a 7 o) = (p Vp) f{uj, \k\) + weitereTerme, 

wobei f(uj, \k\) eine rationale Funktion von uj und |fc| ist. Zwecks Berechnung von P wird dann 
gemafi 

(pVp) f(w, \k\) -> (ip\ (pVp) |V7 f(w, \k\) 
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ersetzt. Die Auswertung der Matrixelemente wird von in Mathematica geschriebenen symbolischen 
Programmen erledigt (vom Autor entwickelt). Bis zur Ordnung a 6 gilt in unserem Fall zum Beispiel 

M (pVp) |V) = a 4 to 3 - — a 6 to 3 . 

Die Auswertung der Matrixelemente ist auch von Hand mit Hilfe von Kommutatorbcziehungcn 
untcr intcnsivcr Ausnutzung der Algebra der Dirac-Glcichung mit Coulombfeld moglich. Dieser 
Weg wurde zur Verifikation der Ergcbnisse cbenfalls bcschritten. 

5. Zusammenfassung des Ergebnisses. Fiir alle r Q 's, deren Koeffizienten nicht Null sind, wird das oben 
beschriebene Verfahren angewandt. Das Ergebnis 

^ = EI>w%miv<) 

0=0 j=0 

wird aufaddiert, vereinfacht und in Potenzen der Feinstrukturkonstante a geordnet. 

Den letzten Schritt bei der Auswertung des Hochenergie-Anteils bilden die Integrationen beziiglich |fc| 
und u>. Die k- Integration kann mit Residuenintegration ausgefuhrt werden. 

Fiir die oj-Intcgration werden ultraviolctt divergente Ausdriickc regularisiert, es werden Feynman- 
Parameter eingefiihrt und die Integration kovariant ausgefuhrt. 

Fiir die ultraviolett konvergenten Ausdriicke wird zunachst die Intcgrationsvariable rcskaliert und auf 
einc dimensionslose Grofie zuruckgefiihrt. Bei der cj-Integration beschrankt man sich auf den oberen 
Ast des C#-Weges (das Gesamtergebnis ist gerade das Integral iiber den oberen Ast, plus das komplex 
Konjugierte dieses Integrals). 



3.8 Der Niedrigenergie-Anteil 



Bei der Behandlung des Niedrigenergie-Anteils wird die w-Integration als erstes ausgefuhrt. Wir gehen 
aus vom allgemeinen Ausdruck fiir den Niedrigenergieanteil zur Selbstenergie 

J Cl (2tt) 4 p- jk — m- TV 

Zunachst formcn wir das Matrixelement mit dem Propagator in einer fiir den Niedrigenergieanteil 
giinstigen Weise um. Untcr Ausnutzung der Bezichung 

1 



A±C = 
B C^BA- 1 



(3.104) 



und der Identitat 

e ik-f p i e -ik-f = p i + ^ . = p i _ ^ 

die man mit der Campbcll-Bakcr-Hausdorff-Formcl oder durch cxplizitcs Differenzieren nachweist, ergibt 
sich mit der Definition a = 7°7° = I (diesc Definition crfolgt hier rein aus Zweckmafiigkeit und ist nicht 
"wortlich" , im Sinnc cines Lorentz-Index. zu verstchen), 

1 



P, 



OT/ 7 



to — 7°y 
i 



(3.105) 



7°ti; — 7°V" — m)e 



— ik- 



1.4- 1 it ik-r 



7°((i?.0 — u>) — 7 • p — 7°U — m) 



(E^p — lu) — a ■ p— V — 7°m 



a e 



1 

{Eij> — uj) — Hu 



—ik-'. 
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Wir verwenden fur den Niedrigenergieanteil die Coulomb-Eichung. In dieser Eichung ist wie bereits in 
Abschnitt 3.5 erwahnt 



D 0i = Ao = 0, 

Zunachst wollen wir denjenigen Anteil behandcln, der vom Z?oo-Element herriihrt. Es ist 
Mit 

(ri1 (^-o,)-^ |r2> =5 (£-a,)-^(l-i„) (3 ' 107) 
konnen wir den von Dqq hcrvorgerufenen Tcilbeitrag zur Sclbstencrgic wie folgt hinschrcibcn: 

e P • 2 f d3k f dLJ 1 / ,| \*Pn)(lpn\ , 

J {^r Jc L 27r k 2 (Eip - us) - E n (l - it]) 

Da aufgrund der speziellen Wahl des C-Integrationsweges nur diejenigen Pole des Dirac-Coulomb- 
Propagators mit E n < entscheidend sind (der Pol mit E„ = E^ > tragt somit nicht bei), ist 
der Integrand innerhalb des Ci-Weges analytisch. Damit ergeben sich aber bei der w-Integration auf Cl 
(mit uj < e) keine Residuen, die zur Lamb-Shift beitragen konnten. Geht man mit e — ► 0, so verschwindet 
der von der von SElao erzeugte Beitrag. 

Wir schreiben daher 8 El folgendermafien unter Vernachlassigung von Dqq auf: 

J (2tt) 3 J Cl 2ir W 2 _ k 2 v > 

Dabci ergibt sich P^ (u>) unter Ausnutzung von 3.105 wie folgt: 

pv(us) = Wa' ^"j^—^Z —) a3 e-^IV). (3-108) 

P^(us) stellt eine Funktion dar, die innerhalb der Ci-Kontour analytisch ist. Der E — w-Ausdruck im 
Nenner des Elektronpropagators steuert kein Residuum bei, da us < e nach oben beschrankt ist (mit 
w<msi£). Daher konnen wir die (^-Integration mit dem Rcsiduensatz ausfiihren. Es ist 

1 1 



uj 2 — k 2 + ir\ (us — |fc| + iT])(u! + \k\ — irf) 
Der Pol bei us = \k\ — irj tragt zum Integral bei, das Residuum ist 

/d 3 k 1 1 
- 7^3^(- 27ri )^n 5T,iJ ' pii (l fc D ( 3 - 109 ) 
,k\<e ( 2 ^) 2tt 2|ft| 



= — e 



= -e 2 



— S T ' ij a l e lk - r — ± r a'e- lk - r \ip) 



\%\<e (27r) 3 2|fc| H d -(E 4 ,-uj) 
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mit 

u=\k\. (3.110) 
Fiir den Niedrigenergie-Anteil wird daher lu mit |fc| identifiziert. 

Jetzt wird im rclevanten Matrixelement P y eine unitare Foldy-Wouthuysen- Transformation (cigcntlich 
zwei hintereinandergeschaltete unitare Transformationen U und U', siehe Abschnitt 3.4) eingefiigt. Die 
Produkt- Transformation £/fw = UU' wird hier dcr Einfachhcit halber wiedcr mit U bczeichnct. Es gilt 

P« = (V'+kV^- l - -c?e^*\1>) (3.111) 

H D - (E - lu) 

= (^+\U + Ua l e ltF U + U — j- U+ U a 3 e~^ p U + U\ib) 

H D - (E — lo) 

= < W (Ucte*+U+) u{Hd _* e _„ ))u+ (Ua^-U-) \U*) 

= (u<fe**U+) u{Hd _I e _ u))u+ ) (Ua>e-^U+) |fl , 

mit \ip) = U\il)) (bei dieser Rechnung wurde wieder 3.104 benutzt). 

Die Foldy-Wouthuysen- Transformation U hat zum Ziel, die ungeraden Anteile des Hamiltonians bis 
zur Ordnung a 4 , d.h. bis zur crstcn rclativistischcn Korrcktur zum Schrddinger-Hamiltonian, zum 
Verschwinden zu bringen. Wir wollcn noch cinmal das Ergcbnis der FW- Transformation aus dem 
Abschnitt 3.4 (Gl. 3.46) in Erinnerung rufen: 

H FW = UH D U + =U (a-p + pm + V)U- 1 = m + H s + SH (3.112) 

Hpw soil den Foldy-Wouthuysen-(FW)-Hamiltonian bezeichnen. Hierbei ist: 

-,2 
V 

H S = ^- + V 
2 m 

und 

Der Ncnncr in P y kann folgendermafien umgeformt werden (Es = Schrodingcr-Energic, 8E = a 4 - 
Korrektur zur Schrodinger-Energie) : 

U (H D - (E-uj)) U + = H fw -(E-lu) = (3.113) 

= (m + H s + SH) - (m + E s + 5E - lu) - 
= (H s -(E s -oj)) + SH-5E. 

Weil 

Hs-iE s - U ) =G ^ E8 - u) (3 ' 114) 

gcradc dem nichtrclativistischcn Propagator G nr cntspricht. sind wir in der Lage, mit dem (viel 
einfachercn) nichtrelativistischen Propagator zu rechnen und alle relativistischen Korrekturen dazu 
storungstheoretisch zu behandeln. 

Dazu ist es notig, auch die erste relativistische Korrektur zur Schrodinger-Eigenfunktion <p mitzunehmen. 
Selbstverstandlich ist die urspriinglichc Eigenfunktion \<fr) des Schrodinger-Hamiltonians nicht mehr 
Eigenfunktion des FW-Hamiltonians. Mit der ersten relativistischen Korrektur ergibt sich aus Gl. 3.54 

M = |0) + |^) = |0) + — - l _ 5H\4>). (3.115) 
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Das Ergebnis der FW- Transformation der Ausdriicke y 3 = a? exp(ifc • r) ist (siehe Gl. 3.72): 



j 

Vfw 



EL 

m 



1 + i ( k -f 



1 



2 m 3 
2m 



7 -2 



k ■ 



2 m z r" 5 



fc X 



2m 



k x a 



5y 3 



(3.116) 



Wir konnen daher 

Ua j e &f! U + = — 
m 

schreiben mit p 3 /m als Hauptbcitrag und fo/- 7 als den relativistischcn Korrekturen dazu. Der Ausdruck 
P lJ setzt sich somit unter Vernachlassigung hohercr Ordnungen folgcndcrmaficn zusammcn, 

1 



Ua i e ik - p U+ 



U(H D - {E-lu))U+ 



XJa? e 



j e -ik-rjj+ 



m 



(H s - {E s - uj)) + SH -SE 



Z- + 5y 3 

m 



(3.117) 



m Hs — (Es — uj) m 
1 



+2 x (<j>\5y % 



— ^1 

1 



1 



m H s - (E s —w) H s - (E s - uj) m 



1 



1 



m H s - (E s -uj) H s - (E s - uj) m 



+2x (6\E- 



1 



P 



i-\S<jj) + 0(a 4 ) 



m Hs — (Es — uj) m 

Fiir die praktische Auswertung des Niedrigcnergieanteils ist die Einfiihrung einer dimensionsloscn Grofie 
P zweckmafiig, die sich aus P gcmafi 

P = ^S T ' ij P 13 (3.118) 

ergibt. P 13 hat aufgrund des Propagators die Dimension A , m dagegen die Dimension A,^ 1 (in natiir lichen 
Einhciten). Es ergeben sich somit bis zur Ordnung a 2 folgende Beitrage zu P %3 . Zunachst lautet der 
Hauptbeitrag 



Pna 



1 



2 m Hs — (Es — uj) m 



(3.119) 



3m Hs — (Es — uj) 



Die Indizes nd stehen hicrbci fiir den "nichtrclativistischen Dipol" . P n d ist der einzige Beitrag zu P in 
der Ordnung a = 1 und damit der einzige Beitrag zum Niedrigcncrgieanteil der Selbstenergie in der 
Ordnung a(Zct) . Allc anderen Beitrage tragen in der Ordnung a 2 zu P und damit in der Ordnung 
a(Za) 6 zur Lamb-Shift bei. 

Diese weiteren Beitrage wollen wir jetzt der Reihe nach auflisten. Zuerst wird der Beitrag Ps v , der durch 
die FW- Transformation der ct^-Matrizen zustande kommt, aufgcfuhrt, 



p Sy = 2 x %F«{4>W 1 —± , p -\ 

2 Hs - (E s - uj) m 



(3.120) 
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Weitere Beitrage ergeben sich aufgrund der relativistischen Korrektur zum Energieeigenwert, 

Pse = ^S T ' ij {ct>\— ± -5E- 1 -£\4>) (3.121) 

2 m Hs — (Es — w ) Us — (Es — uj) m 

= ^(^77 A 



3m xr '^ H s - (E s -u) H 8 - (E s - u) * ^ n 
aufgrund der relativistischen Korrekturen zum Schrodinger-Hamiltonian, 

Psh = r SH — --— - £\4>) (3.122) 

2 m Us — (Es — <jJ) Hs — (Es — u>) m 

= A \ 5H u A rfW* 

3m Hs - (E s -u) H s - (E s - u) 

und aufgrund der relativistischen Korrektur zur Encrgic-Eigcnfunktion. 

Ps* = 2x^^(01^- 1 (3.123) 

2 m Hs — (Es — to) m 

= ^^ pl H W 

3m Hs — (Es — u>) 

Den letzten Punkt bei der Bchandlung des Niedrigcncrgieantcils bildct die d 3 fc-Integration. Fiir den 
Nicdrigcncrgie-Anteil ist |fc| = oj, daher gilt 

E L = -e 2 / ^ * 6 T >*PV(\k\) (3.124) 
J\k\<e (2tt) 3 2\k\ 

9 f d 3 k I 
J\k\<t ( 27T ) m\k\ 

J\k\<e ( 27r ) m\k\ 



2a 
7rm 

Wir fiihren nun eine Variablcntransformation gcmafi 



(1ujluP(uj). (3.125) 



yf^hnEl _ (am)/2 
^/-2m(E 4> - lu) ~ ^/-2m(E 4> - to) 
zur dimcnsionslosen Grofie t durch. Wir verwenden die Abkiirzung 



(3.126) 



Damit folgt 



mit 



A=^. (3.127) 
dwwP(t(w)) = ^ f dt^-P(t) (3.128) 



A a 
y/X 2 + 2me a 2 + 8e/m' 



(3.129) 



F ergibt sich dann zu 



If 1 l-t 2 

F=-- dt—^P(t). (3.130) 



Wie wir in Abschnitt 4.2 sehen werden, ist der Ubergang zur dimensionslosen Variablen t auch dadurch gut 
motiviert, dafi die Propagatoren in der gebrauchlichcn Darstcllung die Energie gerade mit einer zu t direkt 
proportionalen Variablen parametrisieren. t ist also in gewissem Sinne eine "natiirliche" Parametrisierung 
von E. 
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Kapitel 4 

Die Selbstenergie des 2P-Zustands 
(j=l/2) zur sechsten Ordnung in 
Z_alpha 



4.1 Der Hochenergie-Anteil 
4.1.1 Allgemeines zur Rechenmethode 

Bereits im Abschnitt 3.7 war gezeigt worden, da8 man den Hochenergie-Anteil in Potenzen des bindenden 
Coulomb-Potentials aufspalten kann. Bis zur Ordnung a e tragen zum Matrix-Element 

P = ^ u v I ov^M 

p— f — m — 7 U V 

gerade der 0-,l-,2- und 3- Vertex- Anteil bei. Wir wollen uns noch einmal die Form der Matrix-Elementc 
am Beispiel des 2- Vertex- Anteils in Erinnerung rufen. 

P2 = (^|7V^7 ^7°V^7^) (4.1) 
= $|M 2 |V>) 

mit der Matrix 

M 2 = 7 r^l°V^ V^. (4.2) 

Hier soil auf die rechcntcchnischen Details der in Kapitel 3.7 beschriebenen Methode zur Auswcrtung 
eines Vertex- Anteils nalier eingegangen werden. 

1. Entwicklung der Mj-Matrizen in die 16 T-Matrizen gemafi der Gl. 3.24: 

P 

Wclche Matrizen hier beitragen. hangt vom Matrix-Element ab. Fur den 3- Vertex- Anteil bekommen 
wir Beitrage nur von der Id-Matrix und von der 7°-Matrix, fur den 2- Vertex tragen die Id-Matrix, 
die 7°-Matrix, die c^-Matrizen und auch die E-Matrizen bei. 

2. Entwicklung der a/3-Koeffizienten in Potenzen der Feinstrukturkonstanten a. Diese Entwicklung 
bcruht im wescntlichcn auf folgcndcr Bcobachtung: Die Matrizen sind aus frcicn Fcynman- 
Propagatoren 

F i>— }( — m 
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mit dazwischen liegenden Potentialtermen aufgebaut. Den freien Feynman-Propagator kann man 
auch schreiben als 

_ i>— ft + m 
(p — fc) 2 — m 2 

Mit dcr Bindungscncrgie E b des Elektrons 

E b = E S + 5E + ... 

(Es = Schrodinger-Energie und SE = crstc rclativistische Korrektur) lafit sich der Nenner umformcn 
zu 

(p — fc) 2 — m 2 = (m + Eb — oj) 2 — (p— fc) 2 — m 2 

= m 2 + 2mE b + E 2 + 2mu + 2E b uj + lu 2 

—p 2 + 2k ■ p — k 2 — m 2 

= (lu 2 - k 2 - 2muu) + (2mE b + E 2 + 2E b uo - p 2 + 2k ■ p). 
Wir erkennen an diesem Ausdruck den a°-Haupttcil 

lo 2 — k 2 — 2mu> 

sowie hohere Terme 

2mE b + E b + 2E b cu - p 2 + 2k ■ p. 

In diesen hoheren Termen wird nun bis zur notwendigen Potenz in a entwickelt. Dadurch "rutschen" 
die Impulsopcratoren des Elektrons p und Ausdriicke mit den Photoncn-Impulsen k-p in den Zahler, 
und es entstchen neuc Matrix-Elemcntc, die auszuwerten sind. Bei der Entwicklung mufi man 
etwas Vorsicht walten lassen, denn die verschiedenen Impuls-Operatoren, die von den einzclncn 
Feynman-Propagatorcn kommcn. sind zu untcrschcidcn, da die von ihncn crzeugten Matrix- 
Elemcntc untcrschicdlichc Wertc haben. 

3. Mittelung iiber die raumlichen Anteile des Photonenimpulscs. 

Wie wir im vorigen Abschnitt geschen haben, kommcn bei dcr Entwicklung des Nenners des 
Elektroncnpropagators Terme der Gestalt k ■ p vor. Da wir jedoch spater iiber den gesamten R 3 
integricrcn wollcn, konncn wir an diescr Stcllc bercits die Mittelung iiber die raumlichen Winkcl 
vornehmen und dadurch die Ausdriicke stark vereinfachen. 

Die Mittelung iiber die Winkcl crfolgt durch die Ersctzungs-Vorschrift: 

- ^ / dnf(k). (4.4) 
Fiir diese Mittelung werden unter anderem die folgenden Formeln benotigt: 

1/3 k 2 (p-q) (4.5) 

l/h(k 2 ) 2 p 2 (p-q) (4.6) 

1/15 (fc 2 ) 2 (2(p.?f+pV) (4.7) 

1/7 (Pf(p-pf (p-q) (4.8) 

1/35 (fc 2 )' (p- p) (A(p- q) 2 +p 2 q 2 ) (4.9) 

1/35 (fc 2 ) 3 (2(p-q) 2 +3p 2 <f). (4.10) 



-J d n(k-p)(k-q) = 

±-JdCl(k-p) 3 (k-q) = 

i- J dil(k-p) 2 (k-q) 2 = 

i- J dn(k-p) 5 (k-q) = 

i- J dil(k-p) i (k-q) 2 = 

±-J d n(k-p) 3 (k-q) 3 = 
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Diese Formeln gelten fur zwei beliebige Vektoren p und q. 
4. Ersctzung dcr Matrixelemente. In Abschnitt 3.7 war die Ersctzung dcr Matrix-Elcmente am Bcispicl 

(pVp) f(w, \k\) -> (V| (pVp) |V) /(w, |fc|). 
crlautcrt wordcn. Bis zur Ordnung a 6 ist das gesuchte Matrix-Element gegeben durch 

(VI (pVp) lib) = — — a 4 m 3 - -^-a 6 m 3 . 

Zur computergestiitzten Berechnung dieses Matrix-Elements wird die Wellcnfunktion des 2P\/2- 
Zustandes explizit cingegeben und bis zur Ordnung a 6 entwickclt. Die Impuls- und andcre 
Opcratoren werden im Ortsraum explizit auf die Wellcnfunktionen angewandt. Das Rcsultat 
wird in Potenzen von a geordnet. Dabci cntstchen Ausdriicke mit mchrcren tausend Tcrmen. 
Alle diese Terme haben allerdings eine Standard-Struktur. so daB die Integration (bis auf einige 
Ausnahmcn) mit wcnigcn Intcgrationsregcln durchgcfiihrt wcrdcn kann. Dadurch wird die Rcchnung 
um Grofienordnungen beschlcunigt. Wiirde man die in Mathematica eingebauten Integrations- 
Routinen benutzen, so wiirden die Rcchnungen auch auf den schnellsten Maschinen bei weitem zu 
lange dauern. 

Die Auswertung der Matrixelemente ist auch von Hand mit Hilfe von Kommutatorbcziehungcn 
untcr intensiver Ausnutzung der Algebra der Dirac-Gleichung mit Coulombfeld moglich. Bei Matrix- 
Elementen, die nur in der a 6 -Potenz beitragen, ist iibrigens die Verwendung der nichtrelativistischen 
Approximation an die Dirac-Wellenfunktion \ifj), das heiBt der Schrodinger- Wellcnfunktion \<j>), 
moglich. Hicr wollen wir nur einige sehr haufig verwendete Beziehungen unter den Operatoren 
vorstellen, die zur Auswertung der Matrix-Elcmente von Hand niitzlich sind. Es sind dies die fur 
IV) geltendc Dirac-Gleichung 

a ■ p '= E — 7°m — V, 
die fiir \<j>) geltende Schrodinger-Gleichung 

p 2 = 2m(E - V) 

sowie die Kommutatorbezichung 

ABA = 1/2{[A, [B, A]] + A 2 B + BA 2 ). 

Diese Kommutatorbeziehung kann man zum Beispicl anwenden, um das Matrix-Element 
("01 ipVp) folgendermafien umzuformen: 



(V| (pVp) |V> - (VI (-A(V) + 2fv) |V). 



Dabei ergibt A(V) einen der 5-Funktion proportionalen Term. In erster Ordnung tragt dieser Term 
nicht bei, da die P- Wellcnfunktion am Ursprung vcrschwindct. Die erste relativistischc Korrcktur 
indes verschwindet nicht, so daB sich in Ordnung a 6 ein weiterer Beitrag ergibt. 

5. Zusammenfassung des Ergebnisses: Nach der Ersetzung aller Matrix-Elcmente wird das Ergcbnis 
in Potenzen von a zusammengefafit und geordnet. Alle Vertex- Antcilc tragen zur Ordnung a 6 bei, 
da wir einen gebundenen Zustand bchandcln. 

Zuletzt eine Bemerkung zu den angegebenen Teilergebnissen: In alien Zwischenergebnissen fur die Matrix- 
Elemente und auch fur die Resultate zu den einzelnen Vertex- Termen ersetzen wir der Einfachheit halber 
den Faktor Za durch a, um ihn dann im Endergebnis voll auszuschreiben. 



4.1.2 Der 3- Vertex Anteil 

Der 3- Vertex- Anteil (3 Vertizes des Coulomb-Feldes zwischen den Photonen-Vertizes) wird durch das 
Feynman-Diagramm 4.1 wiedergegeben. Da dieser Beitrag am einfachsten auszuwerten ist, wollen wir 
uns ihm zucrst zuwenden. Es ist fiir den 3- Vertex- Anteil das Matrixclcmcnt 
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Abbildung 4.1: 3- Vertex- Anteil 
^3 = (Mfl^fV^V^V^^) (4.11) 
= &\MM 

Ah = 7V^7V^ 7 V^ 7 V^7m (4-12) 

auszuwerten. Zunachst zur Entwicklung von M3 in die 16 T-Matrizen. Bis zur Ordnung a 6 sind nur die 
Entwicklungs-Koemzienten von Id und 7 von Null verschieden. Der Entwicklungskoeffizient von Id lautet 

a w = V 3 [2(-(k 2 ) 2 -8k 2 m 2 -8m i + 12k 2 mu + 16m 3 uj- (4.13) 



mit der Matrix 



-6k 2 u 2 - -\2m 2 uj 2 + 4mu 3 - w 4 )^j / (k 2 + 2 m lj - uj 2 ^ 



Fiir den Entwicklungskoeffizienten von 7 ergibt sich 



16 (fc 2 m 2 +2m 4 -fmw-4m 3 w + 3ra 2 w 
V 3 i 



mcu 3 



(P + 2muj - uj 2 ^ 



Sowohl aid wie auch a 7 o enthalten aufgrund des V 3 -Terms ausschlicBlich Antcilc in a 6 . Da keine 
Entwicklung des Nenners des Elcktronenpropagators in den Kocffizienten notig ist, enthalten die 
Kocffizicntcn keine Skalarprodukte von k mit anderen Vektoren, und wir brauchen keine Mittclung iibcr 
die Winkel vorzunehmen. Die Koeffizienten liegen bereits in der Form 

V 3 f(u;,\k\) 

vor. Wir konnen also die Ersctzung bezuglich der Matrixclcmente gemafi 

a ld = V 3 f(u>,\k\) ^ {iP\V 3 \i>)f(u>,\k\) 

und 

a 7 o = V 3 f{u, \k\) -> (iP\ 7 V 3 \iP)f(u, \k\) 
vornehmen. Die Matrixelemente lauten cxplizit (bis zur Ordnung a 6 ): 

Wv 3 W = Wi°v 3 W = -^« 6 ™ 3 - 

Die besondere Form der P-Zustand-Wellenfunktion bedingt, dafi die obigen Matrixelemente in a 6 
konvergieren und somit die Auswertung durch eine Entwicklung des Ausdrucks 3.88 in Potenzen des 
Coulomb- Feldes m5glich ist. Damit ist das Endergebnis fiir den 3- Vertex- Anteil 

P 3 = a 6 (^{k 2 ) 2 m 3 -8m 7 -4:k 2 m i u> + 16m 6 uj + 6k 2 m 3 Lu 2 - (4.14) 



-12m 5 w 2 +4m 4 w 3 +m 3 uj 4 ^ j f 12 [k 2 + 2 m u - w 2 ) 



Somit ist der einfachste der Vertex- Anteile berechnet. Es gibt fiir den 3- Vertex- Anteil nur einen a 6 - 
Kocffizicntcn, da hicr bereits 3 Coulomb- Vertizes enthalten sind, die einen Faktor (Za) 6 bcitra gen. 
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4.1.3 Der 2- Vertex Anteil 

Der 2- Vertex- Anteil ist durch das Feynman-Diagramm 4.2 gegeben. Zwischen den Vertizes der virtuellen 




X X 

Abbildung 4.2: 2-Vertex-Anteil 
Photonen wirkt das Coulomb-Feld zweimal auf das Elcktron. Es ist somit das Matrixclcmcnt 

P 2 = M/yl/ylyyl^) (4 . 15 ) 

= {4>\M 2 \i>) 

mit der Matrix 

M 2 = 7 V-l 7 0yl/^I 7M (4.16) 

auszuwerten. Bei der Entwicklung von M 2 in die T-Matrizen erhalten wir weit grofiere Terme als beim 
3- Vertex- Anteil. Dafiir sind zwei Griinde verantwortlich. Zum einen erhalten wir nichtvcrschwindendc 
Kocffizicntcn sowohl fur die Id-Matrix, die 7°-Matrix, die a'-Matrizcn wic auch fiir die S-Matrizcn 
(allc andcrcn Projektionen vcrschwindcn in der Ordnung (Za) 6 ). Zum anderen muB der Nenner im 
Elektroncnpropagator in a cntwickelt wcrdcn, was zu grofiercn Termen fiihrt. 

Aufgrund der Lange der auftretenden Terme verzichten wir auf eine Darstellung der Rcchnung und gehen 
gleich zu den Matrix-Elementen iiber. Fiir die Ersetzung durch Matrix-Elemente in den KoefRzienten 
werden folgende Ausdriicke benotigt: 

//iTr2in a 4 m 2 47a 6 ra 2 f .+~\ 



a 6 m 4 



(VWV) = -^g-, (4.18) 



6 4 

a m 



{^\VpVm = (4.19) 

ty\V? 2 V\il>) = (4.20) 

c 6 4 
i Of 111 

(i>\?V 2 m = —g-. (4.21) 

Zur Auswertung der Matrix-Elemente werden im Abschnitt 4.1.1 einige Erklarungen gegeben. Fiir den 
7°-Beitrag benotigcn wir das folgende Matrix-Element 

<#y°V fl |*> = ^ + ^. (4.22) 

Alle anderen Matrix-Elemente fur den 7°-Koemzientcn stimmen bis zur Ordnung a 6 mit den 
entsprechenden Matrix-Elementen fiir den KoefRzienten von Id iiberein. Fiir den KoefRzienten der a- 
Matrizen benStigen wir folgende Matrixelemente, 

{i>\V 2 {a-m) = -g— , (4-23) 

^\V{a-p)V\^) = ^— . (4.24) 
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SchlieBlich benotigen wir fur den Koeffizienten der E-Matrizen noch folgende Matrix-Elemente: 



6 ™4 



(iP\V£-(pV)xp\iP) = —a 



(i/j\Y, ■ (pV 2 ) x p\tp) = — a e m 4 
Das Endresultat fiir den 2- Vertex- Anteil lautct: 

P 2 = a 4 ( k 2 m 3 - 4 m 5 - 3 k 2 m 2 w + 6 m 4 u - 



(4.25) 
(4.26) 



3raW-ra^w d / 6 fc +2mw-u/ + 



a b [12.2 {ky m 3 + 492 {k z y m b - 432 r m' - 141 (fc^) 3 w 

-340 (fc 2 ) 2 m 4 w + 968 k 2 m e lu + 672 to 8 cj - 558 (fc 2 ) 2 to 3 cj 2 - 

-264 fc 2 m 5 w 2 - 1904 to 7 lu 2 + 235 (A? 2 ) 2 m 2 uj 3 - 

-644 fc 2 m 4 uj 3 + 2120 to 6 w 3 + 280 fc 2 m 3 uj 4 - 1052 to 5 cj 4 - 



47k 2 m 2 to 5 + 88m 4 w° + 146m J w D - 47m z w 



,4, ,5 



3 , ,6 



144 Ik +2muj - ui' 



(4.27) 



4.1.4 Der 1-Vertex Anteil 



Dem 1-Vertex- Anteil (1 Vertex des Coulomb-Feldes zwischen den Photonen-Vertizes) entspricht dem 
Feynman-Diagramm 4.3. Zwischen den Photonen-Vertizes wirkt das bindendc Coulomb-Feld cinmal auf 




Abbildung 4.3: 1-Vertex-Anteil 
das Elcktron. Es ist somit das Matrixclcment 



(4.28) 



auszuwerten, wobci fiir die Matrix M\ gilt: 



(4.29) 
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Wie beim 2- Vertex- Anteil verzichten wir wegen der Lange der Termc auf eine Auffiihrung der 
Entwicklungskoeffizienten. In der Ordnung a e erhalten wir Beitrage fur die Id-Matrix, die 7°-Matrix, 
die a l -Matrizen wie audi fiir die S-Matrizen (alle anderen Projektionen verschwinden in der Ordnung 
a 6 ). Anbei eine Liste der fiir diesen 1-Vertex Anteil berechneten Matrixelemente: 



(4>\V\4 

(^|T/(p 2 ) 2 |^ 
(ip\p 2 Vp 2 \tp 

(ip\ (pVp)p 2 \i> 



-a 2 m hoftm 63 a 6 m 



4 32 512 

-5a 4 m 3 337 a 6 m 3 



48 



1152 



48 

-35 a 6 m 5 
192 1 

-19 a 6 m 5 
192 : 

-9a 6 m 5 
64 ' 

—5 a e m 5 



1152 



64. 



(ip\ (p(pVp)p) \i 

Fiir den 7°-Koeffizicnten benotigen wir zusatzlich weitere Matrix-Elemente 

(^ItVIV) 



(^Vp 2 \i>) = 

(4>\f{pVp)p 2 \4>) = 
(iPh°(p (pV0)p)\i>) = 



4 

—5 a 4 m 3 
48 

—5 a m 
48 

-35 a 6 m 5 
192 : 

— 19 a 6 m 5 
192 : 

—9 a e m 5 
64 ' 

—5 a 6 m 5 
64 ' 



16 32 

43 a 6 m 3 
288 ' 

113 a 6 m 3 
576 ' 
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Fiir den Kocffizicnten der a l -Matrizen erhalten wir die folgenden Matrixelemcntc: 

—7 a 4 to 2 65 a 6 m 2 



(fl>\V(3-0)\1>) 
■p\V(a-p)V\TP) 



48 288 ' 

-37 a 6 m 4 



192 



23 a 6 m 4 



(a -p) (p Vp) |V) = 



192 

Schlicfilich bcnotigen wir fiir den Kocffizicnten der S-Matrizcn noch folgende Matrixelemcntc 

i 1 91 ?' 

<V>|£-(i?F)xi#) = -a 4 m 3 + — Vm 3 , 

12 57o 



a 6 m 5 . 
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Das Ergebnis fur den 1- Vertex- Anteil lautet: 

a 2 (k 2 m — 2m 3 + 2m 2 uj + muj 2 ) 

Pi = —± - 2 L + (4.50) 

2[k 2 + 2muj -uj 2 \ 

+a 4 [45 (fc 2 ) 3 m + 46 (fc 2 ) 2 m 3 - 48 k 2 to 5 + 188 (fc 2 ) 2 m 2 w - 



-256 k 2 m 4 oj + 384 to 6 uj - 45 (A: 2 ) 2 to uj 2 + 
+576 fc 2 m 3 uj 2 - 888 m 5 w 2 - 56 fc 2 to 2 w 3 + 720 to 4 uj 3 - 

-45 fc 2 to w 4 - 102 to 3 w 4 - 132 to 2 uj 5 + 45 to uj 6 ) / f 144 f fc 2 + 2 to w - w 2 



a 6 |^8505 (fc 2 ) 5 to + 19790 (fc 2 ) 4 to 3 - 43296 (fc 2 ) 3 to 5 - 16128 (fc 2 ) 2 to 7 + 

+61420 (fc 2 ) 4 to 2 uj - 22104 (fc 2 ) 3 m 4 a; - 477312 (fc 2 ) 2 to 6 w + 
+430080 fc 2 to 8 uj - 25515 (fc 2 ) 4 to w 2 + 279240 (fc 2 ) 3 to 3 uj 2 - 34224 (fc 2 ) 2 m 5 uj 2 
-474240 fc 2 to 7 uj 2 - 368640 7« 9 - 135000 (fc 2 ) 3 m 2 uj 3 + 1011864 (fc 2 ) 2 to 4 uj 3 
-1040480 fc 2 to 6 w 3 + 1413120 m 8 uj 3 + 17010 (fc 2 ) 3 to w 4 - 561180 (fc 2 ) 2 to 3 uj 4 4 
+2020400 fc 2 to 5 to 4 - 2366880 to 7 uj 4 + 36480 (fc 2 ) 2 to 2 uj 5 - 1269960 fc 2 to 4 uj 5 + 
+2243040 to 6 + 17010 (fc 2 ) 2 muj 6 + 205480 fc 2 to 3 w 6 - 1212480 to 5 w 6 + 
+86360 fc 2 to 2 uj 7 + 280200 to 4 uj 7 - 25515 fc 2 to w 8 + 56670 m 3 uj 8 - 

-49260 m 2 uj 9 + 8505 tow 10 



^34560 (k 2 + 2muj- uj 2 } j . (4.51) 



4.1.5 Der 0- Vertex Anteil 

Im 0- Vertex- Anteil rechnen wir mit der Selbstenergiefunktion des freien Elektrons, das heifit, wir werten 
das Diagramm 4.4 aus (die aufieren Spinoren sind wie in jeder Berechmmg der gcbundenen Selbstenergie 
die Dirac-Coulomb-Wcllenfunktioncn, nicht die freien Spinoren). Es ist somit das Matrixelement 



Abbildung 4.4: 0- Vertex- Anteil 
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Po = (V^|7V^7mIV'> (4.52) 

= $\M 3 \1>), (4.53) 

auszuwerten, wobci fiir die Matrix Mq gilt: 

M = 7 V^7/ 1 - (4.54) 

Wie iiblich beschranken wir uns auf die Angabe der relevantcn Matrix-Elcmcnte. Wir crhaltcn Beitragc 
fiir die Id-Matrix, die 7°-Matrix und die a*-Matrizen (alle anderen Projektionen verschwinden in der 
Ordnung a 6 ). Fiir den Beitrag der Id-Matrix ben5tigen wir die folgenden Matrix-Elemcnte 

(ip\ip) = 1 (Normierung), (4.55) 

/ § ,2 1 , , a 2 m 2 13 a 4 m 2 23 a 6 m 2 „. 

WW = — j-+ 48 + 72 ' ( 456 ) 

Ml /. 2 \ 2 ln 7a 4 m 4 167 a 6 m 4 .„ „, 

(^l(/)V> - (4.58) 
Fiir die 7°-Matrix werden die nachfolgenden Elemente gebraucht: 

/ ,1 Ol ,\ t ^ 5Q;4 21 I A rm 



19 a 4 to 2 943 a 6 m 2 
' 96 + 4608 



frW) = ^n C '" + ""^no" ■ ( 4 - 6 °) 



,,, o/.2\ 2 l(1 7a 4 m 4 479a 6 m 4 .„ „.,. 

Fiir die Matrix-Elcmcnte mit den a-Matrizen gilt: 

, , . ,_ . . , , a 2 m 5 a 4 to 63 a 6 to „„. 

(tAI (a -p) |V) = -^ + ^2— + -512 - - (4.62) 

7 a 4 to 3 419 a 6 to 3 

/ n\ 2 21 IT 6 TO 5 

W(d-P) (p 2 ) |V) = — (4.64) 

Auf eine Angabe des Endergebnisses fur den 0- Vertex- Anteil wird hier aufgrund der Lange des Ausdrucks 
verzichtet. Der Ausdruck wiirde mehr als eine Scitc fiillen, ist jedoch von seiner Struktur her dem Ergebnis 
fiir den 1- Vertex- Anteil sehr ahnlich. Er ist von der Gestalt: 

Po = ^ 2(m + ^ + (4.65) 

k 2 + 2 TO OJ — LO 2 

W(...)+a 4 (...)+a 6 (...) 

Den a°-Kocffizicnt kennen wir bereits. Er entspricht gcrade dem Intcgranden der parametrischen 
Darstellung des Massen-Gegenterms, die wir in Abschnitt 2.37 kennengelernt hatten. Somit hebt sich 
dieser (bei Integration iiber |fc| und u> divergente) Ausdruck gerade heraus, und das Integral konvergiert. 
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4.1.6 Der Massen-Gegenterm 



Wie leicht zu crkcnnen, wiirde dcr crstc Term im Matrix-Element Pq auch bei bei Integration bzgl. des 
4-Impulses des virtuellen Photons auch nach Regularisicrung des Photoncnpropagators noch divcrgieren. 

Der Massen-Gegenterm (die Massen-Renormierung) gleicht, wie bereits im Kapitel iiber die Theorie des 
Selbstenergie-Bcitrages, diese Divergenz gearde aus. Wir crhalten mit dcr parametrischen Darstellung des 
Massen-Gegenterms 2.37 



= — le 



<Pk 
(2^ 



1 



1 



fc 2 k 2 - M 2 



2(m 



fc 2 - 2mw 



den Integranden des Massen-Gegenterms als 



2(to + w) 
k 2 — 2m uj 



mit k 2 = to 2 — fc 2 . Wenn man dann noch das Matrix-Element 



(V#) = «>| 7 |V>) = 1- 



5 a 4 



21 a b 



8 



128 1024 



einsetzt, so ergibt sich der Integrand des Massen-Gegenterms, (ip\5m\ip) , als Matrix-Element mit 
geschrieben, zu 



(ij\~f 5m\ip) 



-2 (m + uj) 



a 2 (to + uj) 



k 2 + 2muj-uj 2 A (k 2 + 2rnuj -uj 2 } 

5 a 4 (to + w) 21 a 6 (m + w) 



64 (k 2 + 2 



TO UJ — LO~ 



512 (fc 2 + 2 



mu> — oj' 



(4.66) 



4.1.7 Die k- und omega-Integration 

Den letzten Schritt bei der Bercchnung des Hochcnergie-Anteils bilden die |fc| und die w-Integrationcn. 
Die Ergebnisse fiir den 0-,l-,2- und 3-Vcrtex-Antcil werden aufaddiert (siehc Gl. 3.95) 

und ergeben das Matrixelement P bis zur Ordnung a 6 . Dieses Matrix-Element wird durch Abzug des 
Massen-Gegenterms rcnormicrt und bezuglich \k\ und u intcgriert. 

Den ersten Schritt bildet die |fc|-Integration. Wir definieren die Abkiirzung 



X — y/2 m oj — oj 2 



(4.67) 



um die Ausdriicke in der Rechnung in ihrer Lange zu kontrollieren. Nach gewissen Umformungen des 
Integranden ist die einzige Formel, die zur Integration bzgl. \k\ benotigt wird 



1 

2nd 



1 



1 



fc 2 - uj 2 (A 2 + fc 2 )™ 

i (-i) n a™- 1 



1 



2M^ 2 + A 2 (n — 1)! d(X 2 ) n ^ 1 2X oj 2 + X 2 ' 
wobei die Betragsfunktion fiir komplexes Argument u> wie folgt erklart ist, 

|w| = sgn(3(w)) u, 



(4.68) 
(4.69) 

(4.70) 
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d. h. 



\w\ = +uj falls 9f(w) > 0, (4.71) 

|w| = -u falls < 0. 

Die Formel 4.68 kann man durch Vorziehen der Differentiation nach X 2 vor das Integral und mit 
Anwendung des Residuensatzes leicht nachrechncn. 

Fiir die aj-Intcgration werden ultraviolett divergente Ausdriickc regularisiert, es werden Feynman- 
Parameter cingcfiihrt und die Integration kovariant ausgcfiihrt. Fiir die ultraviolett konvergenten 
Ausdriicke wird zunachst die Integrationsvariable reskaliert und auf eine dimensionslose GroBc 
zuriickgefuhrt. Es wird definiert 

«=J ± £ (4-72) 
X — % uj 



(mit X = \/2moj — uj 2 ). Die ^-Integration wird durch folgende Uberlegung erleichtert: Der Integrand 
hat eine Verzweigung auf der positiven reellen Achse, d.h. der Realteil des Integranden leicht unterhalb 
der reellen Achse ist gerade das Negative des Realteils oberhalb dieser Achse, wie man an der Formel 
4.68 erkennt, und der Imaginarteil behalt gerade sein Vorzeichen bei. Da auf dem Integrationsweg Ch 
unterhalb der rcllcn Achse in negativer Richtung integriert wird, ergibt sich daraus, dafi der Wert 
des Integrals auf dem untcrcn Ast das konjugicrt-Komplexe des Integrals auf dem oberen Ast ist. Die 
Integration nach u und Entwicklung im Infrarot-Rcgulator e erfolgt mit Formeln, die zum Bcispicl in [19] 
zu finden sind. 

4.1.8 Ergebnis des Hochenergie-Anteils 

Nach der Aufuhrung der \k\ und ^-Integration und Ersctzung von a durch Za lautet das Ergebnis fiir 
den F-Faktor des Hochenergie-Anteils 



F H = - 1 - 
H 6 



4177 103 
21600 ~ 180 ^ ' 



a 2(Za) 2 103 2 
( Za > g~ e 180 *• ' ^' ^ ' 



Die Divergenzen in Termen des Abschncideparametcrs e werden sich durch entsprechende Terme im 
Nicdrigcncrgie-Anteil gerade heraushcben. 

4.2 Der Niedrigenergie-Anteil 
4.2.1 Allgemeines zur Rechenmethode 

In dem gesamten Abschnitt iiber den Niedrigenergie-Anteil wird es darum gehen, ein allgemeines Matrix- 
Elemente der Form 

P allg .H = (0|Oip l ^ rP 3 C 2 \4>) (4.74) 

zu berechnen. Dabei ist Es = die Schrodinger-Energie des Zustands |</>) , und wir fiihren eine Variablen- 
Transformtion zum Parameter t mit 

A 

durch, wobei 



A = — sowie X = J-2m(E (f> - uj) (4-75) 

(siehe Abschnitt 3.8, insbesonders Gl. 3.126). 0\ und O2 sind Operatoren, die den Drehimpuls nicht 
andcrn, also beispielswcise skalarc Funktioncn von r. 

Haben wir P(oj) bzw. P(t) erst einmal berechnet, so ergibt sich der gesuchte Beitrag zu F durch 

If 1 1-t 2 
F =-2, d > — P W 
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(siehe Gl. 3.130). Wie ein solches Matrix-Element, welches den Propagator beinhaltet, im einzemen 
berechnet wird, wollen wir hier skizzieren. 

Wir beschranken uns in dieser allgemeinen Diskussion auch auf den Fall der Darstcllung des Propagators 
im Ortsraum und mit Laguerre-Polynomcn. Wie wir im Abschnitt 4.2.2 noch sehen werden, wird 
fur manche Matrix-Elemente eine andere Darstellung des Photonenpropagators gewahlt. Bis auf ein 
einziges Matrix-Element wird jedoch im Ortsraum gerehnet, und wir haben als Ansatz fur den 
Photonenpropagator 

G(r 1 ,r 2 ,E)=J29i(ri,r2,v)Y lm (6 1 ,ct> 1 )Y l * m (62,cj>2) (4.76) 
im 

mit noch zu bestimmendcn Funktionen gi(r\, T2, v). Die Ausdriickc fiir giirx^r^^v) bcinhaltcn noch cine 
Summation iiber Laguerre-Polynome fc-ter Stufe von k = bis oo (siehe Abschnitt 4.2.2). Die Bcrechmmg 
von P a iig. erfolgt nun in folgcndcn Schritten: 

1. Ausfiihren der Winkelintegration. Da die P-Wcllenfunktionen den Drehimpuls / = 1 haben und 
durch die Wirkung von p 1 in Komponenten fur I = und I = 2 aufgespalten werden, ist die 
Winkelintegration bcziiglich der Kugelfunktioncn rclativ einfach durchzufiihrcn. Einzclhcitcn dazu 
finden sich im Abschnitt 4.2.2, in dem auf die mathematischen Grundlagen der Auswertung von 
P-Matrix-Elcmcntcn cingegangen wird. 

2. Ausfiihren der Integration des Radialtcils. Die verblcibcnde Integration beztiglich ri und r-i wird 
mit Formcln erlcdigt, die cbcnfalls in Abschnitt 4.2.2 diskuticrt werden. 

3. Summation iibcr k. Dicsc Summation wird mit Formcln erlcdigt, die zum Toil durch computcrgesttt- 
zte Rechnungen hergeleitet worden sind. Alle so nicht in der Literatur vcrzcichncten und cigens 
hergeleiteten Summcnformeln wurden intensiv an Bcispiclcn gepriift. 

4. Zusammcnhangsformcln zur Vcrcinfachung der hypcrgcometrischen Funktionen. Durch die 
Anwendung von Gaufischcn Zusammcnhangsformcln der hypcrgcometrischen Funktionen werden 
diese vereinfacht, und das Ergebnis wird in eine "lesbare" Form gebracht. 

Durch Integration nach t gemafi Formel 3.130 ergibt sich dann der gesuchtc F-Faktor. Dicsc 
Integration erfolgt in folgcndcn Schritten: 

(a) Die hypcrgcometrischen Funktionen werden um die Stcllc u> = (t = 1) cntwickelt, um eincn 
bei t — 1 konvcrgenten Ausdruck zu crhaltcn. Der cntstchcndc Ausdruck (der "divergente 
Antcil") cnthalt alle Divcrgcnzcn fiir t — > (u> — > oo). Der verblcibcnde ( "nicht-divergente" ) 
Ausdruck ist bei t = 1 und t = konvcrgent und liefcrt cinen weiteren Beitrag zum F-Faktor. 
Da dieser Ausdruck hypcrgcomctrischc Funktionen cnthalt, wird er numcrisch integriert (es 
handelt sich um eine einfachc cindimensionale Integration). 

(b) Integration des divergenten Anteils mit Mathematica-eigenen Routinen von t = t e bis t = 1 
(mit t e gegeben durch 3.129). 

(c) Symbolischc Entwicklung des Resulats der Integration des divergenten Anteils in a bis a 2 , 
danach in e bis e°. 

(d) Berechnung des nicht-divcrgcntcn numcrischen Anteils als einfachc cindimensionale numcrischc 
Integration (GauBsche Quadratur). 

(e) Addition des divergenten und des nicht-divergenten Anteils ergibt F. 

Da die letztgenannten Schrittc zur Berechnung von F mathematisch nicht besonders anspruchsvoll sind, 
werden sie in dieser Arbeit nicht naher erlautert. Wie bereits im Hochenergie-Anteil ersetzen wir in den 
Zwischenergebnissen der Einfachhcit halber den Faktor Za durch a, um ihn dann im Endergebnis voll 
auszuschreiben. 

4.2.2 Mathematische Vorbereitungen 
Verschiedene Darstellungen des Propagators 

Fiir unsere Rechnungen brauchen wir haufig den Propagator 
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der Schrodinger-Gleichung. Im Ortsraum lautet die Gl. 4.2.2 gerade 

(H s - E)G(fx,f 2 ,E) = 5(n - fa). (4.77) 
Man kann diese Gleichung losen mit dem Ansatz (siehe [16]) 

G(n,r 2 ,E) =Y J 9i{rx,r i ,v)Y lm {9 1 Ai)YC m {ez,<j> 2 ) (4.78) 



mit v als "Energie-Parameter" gemafi 



22 1 
v a - 



2mE" 

wobei a = a Bohr = I /(am) der Bohrsche Radius ist. Setzt man nun fur unsere Arbeit 

E = E^ - w, 

so folgt nach einer kurzen Rechnung 

im 

wobei der i-Paramctcr gerade der in Abschnitt 4.2.1 eingcfuhrte Integrationsparameter ist 

Fiir die Funktion gi(r\,r 2 ,E) (den Radialteil des Schrodinger-Coulomb-Propagators) werden in dieser 
Arbeit zwei Darstellungen benutzt. Zunachst wollen wir uns der Darstellung 

gi{n,r 2 ,v) = 2m ( — ) (rir 2 ) 1 exp ( - Vl +? " 2 ) X (4.80) 
\av J \ av J 

x V - T 2l+l ( T 2l+1 ( ^ 

f-^(2l + l + k)\{l + l + k-v) k \av ) k \av 

mit L h a als dem assoziierten Lagucrrc-Polynom der Stufe (a, b) zuwendcn. 

Den Radialteil des Schrodinger-Coulomb-Propagators kann man auch vermittels Whittaker-Funktionen 
darstellen. 

0(11,1*, v)=mav T(l +l ~ "> J- M„, I+1/2 (**) W v , l+1/2 (^) (4.81) 
(2Z + 1)! rir 2 1 \ av J \ av J 

Fiir zwei der Matrix-Elemente benotigen wir noch die Schwingersche Darstellung des Propagators im 
Impulsraum 



( ) f dpp- T ^- (4.82) 

\2sm7TT / J 1 op 

1 - p A 1 



X ■ 



X 2 {p- p'Y + (pa + X 2 )( P ' 2 + A 2 )ii^ 



mit X = -J— 2m(E ( j } — u>) und r = ma/ X = 2t. 
Winkelintegration 

Da die drei P-Wcllenfunktion die Drchimpuls-Quantenzahl I = 1 habcn, transformieren sie sich unter 
Drehungen wie ein Vektor (siehe zum Beispiel [4], Seite 425 ff.). Man kann daher P x -, P y und P z Orbitale 
idcntifiziercn mit 

(f) l (r) =r l 4>(r) (4.83) 
wobei r 1 = x, r 2 = y und r 3 = z (siehe Konventionen) . Es ist 

a% 771 2 amr 

<P{r) = —^e — (4.84) 

4 V2^ 
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eine Funktion, die nur von r abhangt. Die zu berechnenden Matrix-Elemcntc ergeben sich dann aus der 
Mittelung iiber die drei verschiedenen P-Orbitale. 

Nun folgt eine generelle Eigenschaft von Funktionen des Typs pV- 7 f(r). Aus dem Additionsgesetz fur 
Drehimpulse bzw. ganz explizit aus der Formcl 



pV'/(r) = -i^S ij [(3/(r) + rf'(r))} - i 



fifi _ _gij 

3 



WW) 



(4.85) 



erkennt man, dafi die Funktionen p l r J f(r) aufspalten in cinen Anteil mit I = und einen Anteil mit 
1 = 2. Der I = 0- Anteil ist gerade 



-4^'[(3/(r) + r/'(r))] 



und der I = 2- Anteil ergibt sich aus 



(rf(r)) 



AuBerdem haben wir folgendc generelle Eigenschaft der Kugelfunktionen. Fiir zwei Funktionen hi(9,(f>) 
und h,2(0, </>) m h 

hi{r,6,4>) = Y,ci,m{r)Y l , m {9,<j>) (4.86) 
i.in 



mit der / = Z'-Komponente 



h it i = u(0,4>) = 2jc{/, m (r)li/ tm (fl,<ft) 

m 

(keine Summation iiber /') gilt 

J dfii j dn 2 hl(r 1 ,6 1 ,<j> 1 )J2 Y i',rn(0i,<l>i)Y l 7 tm (92,fo)h$(r,9 2 , ( t>^ 



(4.87) 



(4.88) 



dQ, h* l l=l , (ri,0, <f) h 2 ,i=i' (r 2 , 9, 



Daher gilt mit 
und 



/i(r)=Oi0(r) 



(4.89) 
(4.90) 
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die folgende Umformung: 

P allg .H = (^\0 1 p 1 - 1 -p*0 2 \4>) (4.91) 



1,,,,„ , 1 



l] o lP l — — -v l o 2 u ^ 



3 xr ' ^ Hg-(E 8 -u) 
i J dr 1 dr 2 rlrldn 1 dn 2 ((pV /i)(n, X , 0i))* ^ 5z (n, r 2 , w) x 



X^ m (^, ^(0 2 , <fc) ((p l r< f 2 )(r 2 ,6 2 , <j> 2 )) 
1 



dri dr 2 r\ r\ dill ((pV/i)(n, 0, 0)) ;=o <?o(ri, r 2 , i/) ((pV/ 2 )(r 2 , 0, <£)) i=0 

+ ((pV/i)(n,e,0))f =a ^(n.ra.i/) ((pV/ a )(ra,M)), =a 

dri dr 2 r? r 2 2 (3/i(r x ) + ^(n)) ff0 (n, r 3 , (3/ 2 (r 2 ) + r 2 f' 2 (r 2 )) + 

Y (rifi(ri))g 2 (ri,r 2 ,v) (r 2 / 2 (r 2 )) 
Die Winkelintcgration ist somit ausgefiihrt, und es verbleibt die Integration bcziiglich r\ und r 2 . 

T\ und r 2 -Integration 

Wir wollen uns hier nur mit den Matrix-Elementen beschaftigen, die mit der Darstellung des Schrodigner- 
Coulomb Propagators durch Laguerre-Polynome gerechnet sind. Nach der Winkelintcgration sind folgende 
Formeln fur die Berechnung der verbleibenden Integrale beziiglich r± und r 2 hilfreich. Zunachst ist 



dr r a ln(r) exp(-Ar) = A - T(a + 1) (*(a + 1) - ln(A)) (4.92) 
mit der logarithmischen Ablcitung der T-Funktion 

r'Crl 

= (4.93) 
T(x) 

Eine weitere wichtige Formel ist 

dr cxp(-Ar) r' 1 Z£(r) = (4.94) 



A- 1 ^ 2 F 1 (-n, 7 + l, M +l,i)r( 7 + l)r(^ + n+l) 
n\T(n + l) 

Fiir die Ausfuhrung der r± und r 2 -Integrationen bei Verwendung der Whittaker-Funktionen existieren 
ebenfalls Formeln [18]. 

Summation iiber k 

Im Schrodinger-Coulomb-Propagator tritt bei Darstellung im Ortsraum mit Laguerre-Polynomen (siehe 
Gl. 4.80) die Summe iiber den Summationsindcx k auf. Es trcten dabei zum Beispicl folgende Summen 
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auf: 

E? fc = va-o, ( 4 - 95 ) 



k=0 



^k n e = <f(£ fe ,-n,0), (4.96) 



fc=0 



00 C fe 



E^T = -8^1(1,0,0+1,0- (4-97) 
* — ' a + k a 

k=0 

Die Lerchsche $-Funktion ist definiert als 



00 u 
Z 



fe=0 ( a + fc ) 

Die hypergeometrischc Punktion 2-F1 ist crklart als 

2 F l( a,6,c,,) = £ (4.99) 

mit dem Pochhammcr-Symbol 

Ncbcn dicscn rccht cinfachen Summationen treten in den Berechnungen auch komplizicrtcrc Summen suf. 
Wir haben zum Bcispicl die folgcndc Formcl 

E^T(fc + A) e 2 F 1 (-k, b, c, z) = — ±- x T(A) 2-Pi(A, 6, c, --^-) (4.100) 
k=o (1-0 1 5 

Diese und alle anderen Formeln, die mit dcr Theorie der hypergeometrischen Funktionen 
zusammenhangen, finden sich in [17] und [18]. Mit A = 1 folgt aus der Formcl 4.100 

£ e 2 F X (-*, 6, c, z) = 2 11 ^ ^ (4.101) 

fe=0 * 

Mit dcr Beziehung 

k n+1 i k = £^ (fc"£ fc ) (4.102) 

lasscn sich aus der Formel 4.101 auch Formeln fur Summen der Form 

00 

E^"C fc 2-Fi(-fc,6,c,z) (4.103) 
k=a 

aufstcllcn. Bci dcr Anwcndung von Gl. 4.102 auf Gl. 4.101 treten auf der rechten Seite Ablcitungcn nach 
dem Argument dcr hypergeometrischen Funktion auf. Diesc Ablcitungcn lassen sich mit dcr Formcl 

a *i(l + a, b, c, z) = 2 F 1 (a, b, c, z) + *fr (2*1) (a, ft, c, z) 

beseitigen. Durch computcrgcstiitztc Rcchnungcn wurden so Formeln fur die Summen der Gcstalt (4.103) 
fur n = 1 ... 4 hcrgclcitct, da diese Summon in den spatercn Rcchnungcn vorkommcn. Es crfolgt hier nur 
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eine Angabe des Ergebnisses fur n = 4: 



™ 2 F 1 {l,b,c,-^ i ) 15^(2,6, c,-^) 
£ 2 F 1 (-k,b,c,z) = — (-1 + 2 



k=0 



502^(3,6,0,-^) 602^(4,6,0,-^ 



(i-C) 3 (-1 + 4 



24 2 Fi(5,6, c,-^) 



(4.105) 



mit der Eulerschen Konstanten 7 = 0.577216 .... Mit Hilfe dieser Formeln ist es moglich, Summen wie 

00 

J2k n Z k d2(2F 1 )(-k,b,c,z) 

k=0 

auszuwerten. Es wird dabei zunachst die Ableitung nach dem zweiten Parameter der hypergeometrischen 
Funktion durch den Ausdruck fur den entsprechenden Differenzen-Quotienten ersetzt, 

P / j?\r i, », \ 2Fi(-k,b + e,c,z)- 2 Fi(-k,b,c,z) 

d 2 (2-F1) {-k, b, c, z) -> , 

e 

mit einem e, welches wir spatcr wiedcr gegen Null schicken. Nach der obigen Ersetzung lassen sich 
die hergeleiteten Summations- Formeln anwenden, und im Ergcbnis konnen, wie im nachstcn Abschnitt 
besprochen wird, die Parameter der hypergeometrischen Funktion durch Zusammenhangsformcln gecignct 
verandcrt werden. Zulctzt werden Rcgeln fiir die explizitc Form von hypergeometrischen Funktioncn, die 
bereits im Programm Mathematica implementiert sind, verwandt, um die Diffenrentiation (Ubergang 
zum Differentialquotienten) ausfiihren zu konnen. 

Fiir die folgcnden Summen konnte kcinc geschlosscncn Ausdriickc in Termen hypcrgcomctrischc 
Funktioncn gefunden werden: 

Y^^r 1 -9 2 ( 2 F 1 )(-k,b,c,z). 

k=0 

Diese Summen werden numerisch behandelt. Um Singularitaten bei t = 1 und t = auszuschlicBcn, 
summieren wir nicht von k = aufwarts und multiplizicrcn mit t 2 . Wir definicren die "spcziellc Funktion" 

00 (t=±) k 2 
F abc (t) =t 2 E a 2t + k ° 2 (_fc ' 6 ' C 'TTt^' (4 ' 106) 

k=2 

In unserer Rcchnung kommt diesc Funktion in der Gestalt 

k 



F 122 (t) = t 2 Y, r^ixfc * (- fc < 2 > 2 ' j4t) ( 4 - 107 ) 



k=0 

und 



t—l \ 

t+T „ , „ , , , „ „ 2 



F 366 (t) = t 2 3 L 2 t + k d2 {2Fl) { ~ kl 6 ' 6 ' TT7 ) (4 ' 108) 

vor. Man kann zeigen, dafi Fi 22 (t) und F^ft) fur t — ► gegen Null und fiir t — > 1 schncller als (1— i) 2 gegen 
Null gehen. Da Integrale mit diesen Funktionen numerisch behandelt werden und dafiir eine Gausssche 
Quadratur verwandt wird, die bei glattcn Funktioncn ohne Singularitaten schnell konvergiert, ist dieses 
Verhalten sehr hilfreich. 

Weitere Summationen, die bei der Berechnung der P-Matrixelemente auftreten, involvieren die 
logarithmische Ableitung der T-Funktion, 

*(*) = -ii (4.109) 
T(x) 
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Durch zum Teil computergestiitze Rechnungen lassen sich auch fur diese Summationen Formeln angeben. 
Zwei Beispiele seien hier erwahnt: 

f:r^(2 +fe ) = 7 ^r/ ) (4.110) 

fc=0 ^ 

wicdcrum mit 7 als der Eulcrschcn Konstanten und 



ffc 4 r*(Hfc) = 1 6 ^ 13 ^ 2 32 ^ 3 + (4111) 

7 g (1 + (1 + ioe + e 2 ) 
(-1 + 5 

(l-5£ + lU 2 +£ 3 + 16e 4 ) ln(l-Q 

(-1 + C) 5 ? 

Es gibt jcdoch auch fur dicjenigen Summen, wclchc die "f-Funktion beinhalten, bcstimmtc Falle, fur die 
keine Losung in geschlosscncr Form gefunden werdcn konnte. Es sind Funktioncn dcs Typs: 



00 



2k 

1-t 

l+t 



b-2t + k 

k=2 

Dicsc Funktioncn kommen in unserer Rechnung in der Gestalt 

2/,- 



00 f i=*y 

(*) = e E i-M + k ^ 2 + k) (4 - 113) 



k=2 



und 



00 ( l=£\ 

*6s(*) =* 2 E 3„ + 2 t / +fc ^(6 + fc) (4.H4) 



fc=2 



vor. Auch diese Summen werden numerisch bchandclt. Zum SchluB definieren wir noch die folgenden 
"spezielle Funktionen" , die mit der Lerchschen <I>-Funktion zusammenhangen, 



21; 



die in dicscr Arbeit als 



und 



k=2 



00 



1-t 
l+t 



2 k 



$2i(t) = t 2 Y — i '- j (4.116) 

2 A- 



00 



1-t 
l+t 



* 2s(i) = ' 2 S(^Tif ( " 17) 

vorkommen. Auch diese Funktioncn werdcn numerisch bchandclt. 
Zusammenhangsformeln 

Die hypergcometrischen Funktioncn im Ergcbnis 4.105 haben hohe Werte fur den ersten Parameter. 
Daher crschcint es ratsam, diesen Parmeter durch die bekannten Gaussschen Zusammenhangsformeln zu 



51 



verkleinern (siehe zum Beispiel [17]). Auch durch die Anwcndung der letzten Formel aus den Standard- 
Summationen 4.95 treten hypergeometrische Funktionen mit unangenehmen Werten fur die Parameter auf 
(in diesem Fall sind die 6- und c-Paramctcr zu hoch). Um die Parameter geeignet zu verandern, verwenden 
wir einige Zusammenhangsformeln in einer Gestalt, die das Herabsetzen der Parameter ermoglicht. Es 
sind dies: 

• zum Herabsetzten des zweiten und dritten Parameters mit besonderem Augenmerk auf den dritten 
Parameter 

P , , s (-1 + e) 2 F l {a,-l + b 1 -l + c,z) 
2 F 1 (a,b,c,z) = -i + b (4.118) 

jb-c) 2 F 1 {a,-l + b,c,z) 
-1 + 6 

• zum Herabsetzten des zweiten und dritten Parameters mit besonderem Augenmerk auf den zweiten 
Parameter 



tp i h \ (1-c) 2 F 1 (a,-l + b,-l + c,z) fA nn\ 

2 F 1 (a,b,c,z) = — r (4.119) 

(1 + a — c) z 

^ (-1 + c + z — cz) 2 F 1 (a, b, -1 + c, z) 
(1 + a — c) z 

Diese und die vorige Regel miissen bei der Vereinfachung der Ausdriicke so eingesetzt werden, dafi 
nicht einer der Parameter zu schnell sinkt, da dann einc Vereinfachung im Sinne moglichst kleincr 
Wcrte fiir beidc Parameter nicht mehr moglich ist. 

• zum Herabsetzten des ersten Parameters 

tp i h \ (-1 + a-c) 2 Fi(-2 + a,b,c,z) 

2 F 1 (a,b,c,z) = ~ . (4.120) 

(-1 + a) (—1 + z) 

(2 — 2a + c— z + a z — b z) 2 F\ (—1 + a, b, c, z) 
+ (-1 + a) (-1 + z) ' 

• zum Herabsetzten des zweiten Parameters 

2F 1 (a > b,c,z) = + ( _ i + z) (4.121) 

(2-26 + c- z-az + bz) 2 F 1 (a, -1 + b, c, z) 



(-1 + 6) (-1 + z) 



• zum Herabsetzten des dritten Parameters 



w(nh \ (-2 + c) (-1 + c) (-1 + z) 2 J 1 1 ( a ,6,-2 + c,z) 
2 Fi(a,6,c,z) = — — — — — 4.122) 

(1 + a — c) ( — 1 — o + c) z 

(-1 + c) (-2 + c + 3z + az + 6z - 2cz) 2 F X (a, 6, -1 + c, z) 



(1 + a-c) (-1 -b + c) z 



Abkurzungen 



Da die folgendcn Ausdriicke in dieser Arbeit schr haufig vorkommcn, wollcn wir folgendc Abkurzungen 
vereinbarcn: 

s = — (4.123) 

t + l J 
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sowie 

F{t) = 2 F 1 {l 1 -2t,l-2t,s 2 ) (4.124) 

und 

F 2 (t) = 2 F 1 (l,-2t,l-2t,-s). (4.125) 
Bei diesen Funktioncn ist zur Verbesscrung dcr Konvergenz der dritte Parameter groficr als der zweite . 

4.2.3 Der nichtrelativistische Dipol-Anteil 

Der nichtrelativistische Dipol-Anteil ist durch Gleichung 3.119 gegeben: 

Pnd = ^-Wf-s — A — 

6m Hs — {-as — 

Zur Kontrolle wurde die Rcchnung sowohl im Impuls- als auch im Ortsraum ausgefiihrt. Fiir die 
Rechnung im Impulsraum wird die Fourier- Transformation dcr 2P-Wcllcnfunktion bcnotigt, und es wird 
die Schwingcrschc Intcgraldarstcllung des Schrodingcr-Coulomb-Propagators im Impulsraum verwandt. 

Fiir die Rcchnung im Ortsraum wird die "iiblichc" Darstcllung des Propagators im Ortsraum mit 
Lagucrrc-Polynomen verwendet (siehe Abschnitt 4.2.2). Die Summation beziiglich k wird ausgefiihrt, 
und die entstchenden hypcrgcomctrischcn Funktioncn werden mit Hilfc dcr Zusammcnhangsformcln 
vereinfacht. 

AuBcrdcm wird zum dimcnsionslosen Parameter t [siehe Gl. (3.126)] iibergegangen. Das Ergcbnis fiir 
P n d(t) lautet folgcndcrmaficn: 

256 F(t)t 7 (-3 + lit 2 ) 

Pnd(t) = f 5/ - g - L + 4.126 

9(l-tf(l + tf 

2 t 2 (-3 + 6 t + 3 i 2 - 12 < 3 - 29 t 4 - 122 1 5 + 413 i 6 ) 



9 (-1 + tf (1 + t) 3 



Die Integration beziiglich t mit der Formel 3.130 licfert 

2a 2 

F nd = 0.0400223 + — . (4.127) 

4.2.4 Korrekturen zum Strom 
Der nichtrelativistische Quadrupol F nq 

Im nichtrelativistischcn Quadrupol werden diejenigen Korrekturen zum Strom zusammengefafit, die sich 
aus den 



Vi.fw = — {ik-r) + 



in 
und 

ergeben. Die Rechnung erfolgt im Impulsraum (eine vergleichbare Rechnung im Impulsraum fiir S- 
Zustande findet sich zum Beispiel in [19], und die Rcchnung wird daher hier nicht naher beschrieben) . 
Das Ergcbnis fiir P nq (t) lautet: 

64 a 2 F(t) t 5 (25 - 106 1 2 + 29 i 4 ) 
45 (1 - t) 5 (1 + t) s 



a 2 ^15 - 30 1 - 58 1 2 + 146 1 3 + 812 1 4 + 4630 1 5 - 15174 i 6 
-1418t 7 + 4421i 8 ) / (360 (-1 + tf (1 + tf 
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Nach der Integration beziiglich t ergibt sich folgender Ausdruck fur F nq 

2 49a 2 ln(a) 49a 2 ln(e) 



F nq = -1.20115a 2 + ^f^' ( 4129 ) 



Korrektur aufgrund des p l p 2 -Terms 

Die p^-Korrcktur kommt aufgrund des entsprechenden Terms (siehe Gl. 3.68), 

<™ i /,,//2 ' •••• 

der bei der FW- Transformation der reinen a l -Matrix entsteht, zustande. Der Ausdruck fur den P-Tcrm 
lautet: 

1 „i^2 



3m 3 ^ V J Hs-iEs-uj) 1, 

Die Rcchnung erfolgt im Ortsraum. Das Ergebnis fur P p i p 2(t) lautet: 

P„® = 64Q2 1 ( ^ 5 3 (4-131) 
PP 9(-l + i) 3 (l + t) 3 

64 a 2 F(t) t 5 (-1 - 11 t 2 + 36 t 4 ) 
+ 9(l-t) 5 (l + i) 5 

a 2 1 2 (7 - 14 1 - 23 1 2 - 68 1 3 - 103 1 4 - 302 t 5 + 1271 1 6 ) 
18 (-1 + i) 5 (1 + *) 3 

Fur den .F-Faktor wird schlicfilich erhalten: 

9 Act 2 ln(a) 2a 2 ln(e) 
F p , p 2 = 0.618148 a 2 + ^— ^ (4.132) 

Korrektur aufgrund des r x a- Terms 

Dicsc Korrektur kommt durch den Term 

welcher bei der FW- Transformation der reinen o^-Matrix entsteht, zustande. Es ist das Matrix-Element 

la, ^ .A 1 

— — r 1 — f X (7 



3m 2 ^'V2r 3V ' J H s - (E s - u) 

auszuwerten. Mit ciner Rcchnung im Ortsraum erhalt man fiir P rxa (t): 

64 a 2 F(t)t 5 (— 1 + 5 i 2 ) 

Pry. a ; ~A~, 7d 

9(-l + i) 4 (l+t) 4 



(4.134) 



Das Ergebnis fiir den F-Faktor lautet: 



2a 2 t 2 (-l + 2t + 14i J -47t 4 
9(-l + t) 4 (l + tf 



F rxa = 0.346690 a 2 . (4.135) 
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Korrektur aufgrund des k ■ r- Terms 

Dicscr Bcitrag wicderum kommt durch den Term 

vi, FW = --- + ^(k-^)(kxa) j + ..., 

welcher bei der FW- Transformation des Ausdrucks a* (z/c-r^-Matrix entsteht, zustande. Es ist das Matrix- 
Element 

P, r = 2 X %T« m [±(k.r» (i X g )') ^ _ ^ _ u) g W 

auszuwerten. Eine langere Rechnung ergibt, dafi dieser Beitrag gerade dem — i-fachen von F rx<T (t) 
entspricht. Daher gilt 

F k . r = -0.173345 a 2 . (4.137) 

4.2.5 Korrektur zur Energie 

Bei der Korrektur zur Energie ist lediglich cin Matrix-Element zu bcrechnen, 



Pse = ^(0|P 1 77 t4 \ 5E 77 t4 rfM ( 4 - 138 ) 

3m H s - {E s ~uj) H s - (E s - w) 



Fiir den 2P 1 / 2 -Zustand ist 



SE = — — — OL^vn. (4.139) 
128 



Fiir die Berechnung konnen wir die folgende Beziehung ausnutzen: 



2 d 1 



Hs-iEs-u)) duHs-(Es-Lu) 



(4.140) 



Mit der Bcobachtung, dafi sich der Selbstenergie-Beitrag El aus dem P-Term durch Integration gemafi 
Gl. 3.124 



E L = / dwwP(u) 

7TTO Jo 



und der Bcobachtung, dafi der P-Term dirckt proportional zu El ist, lafit sich somit der P-Tcrm vcrmittcls 
eincr particllen Integration bercchncn. Es entsteht ein Randtcrm und cin wcitcrcr Term, der mit den 
iiblichcn Methodcn im Ortsraum berechnet wird. Das Endergebnis fiir den P-Faktor lautet: 

P = 0.0400946 o? - **^L + ^1^- (4-141) 
48 96 

4.2.6 Korrektur zur Wellenfunktion 

Der Beitrag aufgrund der relativistischen Korrektur zur Wellenfunktion 5(f> ist gemafi Gl. 3.123 gegeben 
durch 

P «* = 5^bS? A \P l \^)- (4-142) 

3m Hs — (Es — cj) 

Es ist die Differentialglcichung 3.53, 

(H -{4>\5H\4>)) \<j>) =(E s -H s )\S(f>), 
zu losen. Deren Losung ist durch die zusatzlichcn Bcdingung 

W4>) = o 
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eindeutig bestimmt. Wir erhalten folgenden Ausdruck fur \S(f)): 



7 3 11 7 o 

q/ 2 777,2 —amr as mz r —amr 
a§ ml r ln(amr) .—amr. 

171 

= + (502 + <503 + (504 + 5cf> 5 . (4.143) 

Um die Rechnung bei dcr praktischen Ausfiihrung iibcrsichtlich zu haltcn, bchandcln wir Pg^ fur jeden 
dcr 5 Antcilc von \5<p) gctrcnnt. Es cntstchcn also 5 cinzclnc Tormc, Ps<j>,i bis Ps<j>,5- 

Fs<f>,i bis Ps(pA konnen im Ortsraum berechnet werden. Die jeweiligen P-Terme enthalten die Funktionen 
F(t) und F^{t) und sind in ihrer Struktur den bereits aufgelisteten Ergebnissen fiir die vorhcrigen Tcrmc 
ahnlich. Dahcr vcrzichtcn wir hier auf cine Darstellung dicscr Ergcbnissc und gchcn gleich zu den F- 
Faktoren iiber. Es gilt: 

F s<j>A = 0.0883825 a 2 (4.144) 
F S 4,,2 = -0.0231015 a 2 (4.145) 

F S ,, 3 = -0.608142 a 2 -1^^ + ^^ (4.146) 

F 5M = -0.0261919 a 2 _ ^1 + ^1 (4 . 147) 

18 36 

Aufgrund des Logarithmus im letzten Term, der zu \d(f>) beitragt, wird die Rechnung fur dicscn Term 
sehr viel komplcxer. Es miissen andere Summcnformcln fiir die Summation iiber k bemiiht werden, und 
bereits das Ergebnis fiir den P-Term wird um einiges langer: 



-512a 2 P 122 (t) t 6 
27 (-1 + t) 2 (l + tf 



4096a 2 # 366 (t) (1 - 2t) (-1 + t) t 6 (l + 2t) 
27 (1 + i) 11 

16a 2 F 2 (i) t 5 
3(1- tf (1 + tf 

2 a 2 -/t 2 (-3 + 6 1 + 3 1 2 - 12 1 3 - 29 1 4 - f 22 1 5 + 413 t e 
9 (-1 + i) 5 (1 + tf 



a 2 t 2 p(t) 



(4.148) 



810 (2-t) (-1 + tf (1+t) 13 



128 a 2 1 7 (-3 + 2 1 + 9 1 2 ) ln(^) 16 a 2 F{t) t 5 q{t) 



9 (1 - tf (1 + tf 405 (-1 + tf{l + t) 5 ' 
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Die Funktionen p(t) und q(t) sind gegeben durch 

p(t) = 1710 + 12825 1 + 34920 t 2 + 25920 t 3 - 38160 t 4 + 391876 1 5 (4.149) 

+2892232 i 6 + 6664816 1 7 + 4416444 i 8 - 15374274 i 9 - 44242888 1 10 

-51115488 1 11 - 37744768 t 12 - 6464172 i 13 + 4472472 1 14 

+4920240 1 15 + 1591846 1 16 + 186385 1 17 



sowie 

q(t) = 135 - 7594 1 2 + 2160 7 i 2 + 23251 1 4 - 7920 7 i 4 (4.150) 

2 1 2 £ 
-2160 1 2 ln( ) + 7920 i 4 ln( ). 

v i + r v i + r 

Diejenigen Terme, welche die speziellen Funktionen F\22(t) und F^lt) nicht enthalten, konncn mit 
den iiblichen Methoden nach i integriert werden. Die rcstlichcn Terme, welche die speziellen Funktionen 
enthalten, werden als numerischer Anteil Ns,f,,5 zusammcngefafit und spater mit anderen numerischen 
Termen gemeinsam ausgewertet (siehe Abschnitt 4.2.8). Da es sich dabei um eine einfache numerische 
Integration handelt, kann das Ergebnis mit hoher Genauigkcit angegeben werden. Als Endergebnis fur 
^0,5 erhalten wir: 

9 2 a 2 \n(a) a 2 ln(e) 
Fs*,b = -0.182817a 2 + ^ + N 54> , 5 . (4.151) 

Zahlt man allc Einzclbcitragc zu Fg^ zusammen, so wird 

, 13a 2 ln(a) 13 a 2 ln(e) 
Ft* = -0.751870 a 2 — ^ + — + N s<t> , 5 . (4.152) 

lo ob 

4.2.7 Korrekturen zum Hamiltonian 

Aufgrund der Form der relativistischcn Korrektur zum Hamiltonian, 

(v 2 ) 2 it a a -* 

8 m a 2m z Am z r a 

konncn wir den Beitrag aufgrund von SH (siehe Gl. 3.122), 

Psh = -^-{^ jf 77? s SH u Tw ( 4 ' 153 ) 

3m H s - (E s -w) H s - {E s - lu) 

aufspaltcn in eincn P$, P p i und P^.s-Term, wobci die Namen den charaktcristischcn Elcmcnten in den 
den 3 Summanden in SH cntsprcchen, das heifit, es gilt: 

Ps = -J~(4>\P 1 77 A \ f ^5(r)) 7? A KP l \<t>) (4-154) 



3m vr '^ H s - (E s - uj) V2m 2 v V H s - {E s - to) 



p p* = -^-(4>\P 1 77 i c {-&■) 77 A T^> ( 4 - 155 ) 

3m Hs — {Es~uj) \ 8m J / Hs — (Es~uj) 



PL - S = -^ H S -(E S -.) H s -(E s -u) m (4 ' 156) 

Die Rechnungen zu diesen 3 Beitragen werden durch die Anwesenheit von 2 Propagatoren in den Matrix- 
Elementen erheblich crschwert. Fur Pg nutzt man die Eigenschaft der (5-Distribution aus und geht 
zur Darstellung des Schrodinger-Coulomb-Propagators vermittels der Whittaker-Funktionen fiber. Fiir 
P p 4 und Pl s benutzt man folgenden Trick: Man schreibt die Operatoren in den Beitragen zu SH als 
Kommutatorcn anderer Operatoren mit Hs — (E s —ld). Dadurch ktirzt sich einer der Propagatoren heraus. 
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Beitrag aufgrund des <5(r)-Terms in 8H 

Das Matrix-Element 



1 ... , 1 /. 2 \ 2 1 



Ps = -^-W-s ^ r (^5(r)) 1 -M) (4.157) 

3m Hs — (Es — \2m^ / Hs — (Es — w ) 

wird im Ortsraum unter Vcrwcndung der Darstcllung des Schrodingcr-Coulomb-Propagators mit 
Whittaker-Funktioncn bercchnct. Das Ergebnis in Termen von t ist: 

- (a 2 t 4 (-3 + 4t-8F 2 (t)t + 7t 2 ) 2 ) 

P s {t) = — i -, 3 '- (4.158) 

27(-l + <) 4 (l + t) 4 

Nach der Integration bcziiglich t ergibt sich als Endcrgcbnis fur F$ 

F S = 0.00165032 a? - + £^E> (4.159) 

6 12 

Beitrag aufgrund des p 4 -Terms in SH 

Dieser Beitrag 

P * = -L^ Hs-Ies-u) [ J i£) H s -{E s -u,) pi M 

oder 

Pp4 ~ 24m 4 Hs - (£ s - w) ^ ^ #s - (£ s - w) 
wird unter Zuhilfenahme der zunachst trivialen Bezichung 

f = 2m{{H s - {E s - U ))-V+ (E s - u))) 

berechnet. Diese Beziehung ergibt quadriert: 

(p 2 ) 2 = Am 2 ((H s -(E s -u;)) 2 + 2{E s ~u;)(H s -{E s -Lu)) (4.160) 

-V(H S - (E s - lu)) - (H s - (E s - u))V + {(Eg -to)- V) 2 ). 

Setzt man diesen Ausdruck fiir p 4 in das Matrix-Element ein, so ergeben sich auf naturliche Weise 3 Terme, 
die zum e°-Koemzicntcn beitragen (der durch (Hs — (Es — lu)) 2 erzeugte Beitrag ist zu vernachlassigen, 
da er als Integral J* duiui P(lu) — ► geht fiir e — > 0). Man erhalt 

1. durch den Term 2 (E s - to) (H s - (E s - w)) 

JW = ^2(E S - w)( ^ 1 _^__ yP V) I (4.161) 

2. durch den Term V (H s - (E s - lu)) + (H s - (E s - lu)) V 

P ^ = - 2 ^ l Hs-(E s --) Vpl ^ (4 ' 162) 

3. durch den Term ((E s - lu) - V) 2 

P ^ = - 2 h^ H 3 -lE s -^ Es - "> - V)2 H s -(Es-^) m - (4 ' 163) 
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Der Term P p 4 3 mufi nun noch weiter bchandelt werden. Fur die zwei Propagatoren gilt, da ((Es — uj) — V) 2 
ein Ausdruck der Ordnung a 4 ist, in erster Naherung 

{{E s -u)-Vf- -J- r= (4.164) 



H s -{E s -Lu) Ky ' ' H s - (E s - u) 

1 1 



H s - {E s - u) - {{E s - u) - V) 2 H s -{E s -Lo) 
Wenn wir die temporare Abkiirzung Q = Es — u cinfuhrcn. so cntspricht die Ersetzung 

H s (l, a)-n = H s -(E s ~ u) (4.165) 

- H s - (E s - w) - ((E s - w) - V) 2 = H s (l', a') - fl' 

im Schrodingcr-Hamiltonian und damit im Schrodingcr-Coulomb-Propagator einer Verschiebung der 
Energie um O — > fl' , der Drehimpuls-Quantenzahlcn I — ► I' und des a — > a' im Potential. Damit konnen 
alle diese Korrekturen storungstheoretisch behandelt werden, und es ergeben sich drei Beitrage P p 4 3a , 
Pp 4 ,3b und P p 4 j3c , die respektive fi — * CI', a — > a' und / — > Z' entsprechen. 

Die Beitrage P p * \ und P p 4 2 enthalten nur einen Propagator und sind von ihrer Struktur her einfach. 
Wir nehmen dahcr gleich das Endcrgcbnis fur die P-Faktorcn vorweg: 

F p 4 ! = -0.302463a 2 + (4 . 166) 

9 9 



F p ,, 2 = -0.00661108a 2 + _ ^(f) . (4167) 



Im Matrix-Element P p 4 3a sind zwei Propagatoren vorhanden, der zwischen den Propagatoren stchcndc 
Ausdruck ist jedoch nur eine skalare GroBe (eine Zahl) und kommutiert mit den Propagatoren. Daher 
konnen wir dicscn Bcitrag nach dcm Muster von P$e berechnen (partiellc Integration, Randterm, dann 
Integration des verblcibcnden Beitrags). Das Ergcbnis ist: 

„ ™„™ 2 Ha 2 ln(a) lla 2 ln(e) , d „„. 

P p 4, 3a = 0.559288 a 2 + ^±-L ^-^ (4.168) 

Schwieriger wird die Auswertung der verbleibenden Terme P p 4 ^ und P p 4 i3c . Zunachst zu P p i^b- Bei der 
Summation beziiglich k treten Terme auf, die der Lerchschen Funktion $ entsprechen. Das Ergcbnis fur 
P p 4 3b ist: 

128 a 2 $(s 2 , 2, -2t) t 7 (-3 + 11 t 2 ) 

P P *36 = — *r 4.169 

P ' 9(1-0 (1+0 

32a 2 P(<)< 5 (5-46i 2 + 57i 4 ) 
+ 9(-l + i) 6 (l + i) 6 + 

a 2 t 2 (-5 + 10 1 - 11 1 2 - 52 1 3 + 453 1 4 + 170 1 5 - 821 1 6 ) 



18 (-1 + i) 6 (1 + i) 4 



Die Integration beziiglich t liefert fur P p 4 3fc : 



, 7 a 2 ln(a) 7 a 2 ln(e) 
F p i,3b = -0.0376386 a 2 ^-L + ^-12 (4.170) 

Den verbleibenden Term P p 4 3c berechnet man mit Hilfe einer Variation des a-Faktors im Schrodingcr- 
Coulomb-Propagator. Das Ergcbnis fur den P-Term ist recht komplcx, in seiner Struktur jedoch dem 
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Resultat fur Ps$,5 ahnlich und kann hier in vcrkiirzter Form angegeben werden. 

P+* = 1024Q2A 2 22(t)t6 a (4-171) 
P ' 3 27 (-l + <) 2 (l + <) 6 V ' 

| 8192a 2 F 366 (£) t 6 (-l + 2t) (l + 2t) 

135 

512 a 2 $ 2 i (t) t 6 
+ 27 (-1 + i) 2 (1 +i) 6 

4096 a 2 $ 23 (t) + t 6 (-l + 2t) (1 + 2t) 
135(1 + t) n 

1024 a 2 *2i(i) i 6 
+ 27(-l + t) 2 (1 +t) 6 

8192a 2 ^ 63 (t) (1 -2t) (-1 + (l + 2t) 
135 (1 + t) 11 

+/(7,F 2 (i),F(i),m(t)), 

wobei /(7 F2(t), F(t), ln(£)) eine Funktion dcr aufgefiihrten Grofien ist, deren Struktur den ensprcchendcn 
Termcn aus Ps<f>.5 ahnclt und die aufgrund ihrcr Lange hier nicht weiter aufgcfiihrt wird. Wichtig ist, dafi 
im Ergebnis cinigc spczicllc Funktioncn enthaltcn sind, die cincr getrennten numcrischen Auswcrtung 
zugefiihrt werden (siehe Abschnitt 4.2.8). 

Als Ergebnis fiir F p i 3c erhalten wir 

F pi , 3c = 0.0895706 a 2 + - + A^ 3 , (4.172) 

mit JVp4 3 C als dem noch auszuwertenden numerischen Anteil, der die speziellen Funktionen enthalt. Somit 
crgibt sich als Endergebnis fiir F p i 

F pi = 0.302146 a 2 + - + N p *, 3c . (4.173) 

Beitrag aufgrund des L ■ S- Terms in 5H 

Der Beitrag 

Pr C = - 

3m H s - (E s - lo) V4m 2 r 3 / H s - (E s - w) 
wird mit dcr Definition 

1 d 

D r = -jj-r (4.175) 
r or 



unter Ausnutzung der Kommutatorbczichung 



[H, D r ] = -^3 - 3 (4.176) 



berechnet. Mit dieser Kommutatorbczichung ergeben sich drci Bcitragc, 
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sowie 

a , . , , 1 



Pl-s,2 = =^a (0b* 77 7™ \ D r (<t ■ L) pfy) (4.178) 

H.q — (kit! — U)) V / 



72m 2 ^'^ Hs-(Es-lu) 

mid 



iy 2 



P*|0>. (4-179) 



H s - (E s - u) 



72m y 1 # s -(E s -w 

Aufgrund der Anti-Hcrmitczitat von D r folgt, daB Pls,i und Pls,i gleich sind, und man erhalt als 
Ergebnis fur den P^.^i-Term: 

8a 2 F 2 (t)t 5 

Pl-sa = — ^7 -3 4.180 

9 (-1 + i) 3 (1 + t) 3 

8 a 2 F(t) t 5 (3 - 74 1 2 + 167 t 4 ) 
27(l-<) 5 (l + t) 5 



v2 +4 



t 4 (-l + 2i-88i 2 -98t 3 + 377i 4 ) 
27 (-1 + i) 5 (1 + i) 3 

Nach Integration bezuglich t erhalt man den F-Faktor 

F L . S , = F W = 0.0103947 a 2 + _ ^Mf) . (4 . 181) 

54 108 

Den verbleibendcn F^.g-Term rechnet man mit einer Variation des Drehimpulses im Schrodinger- 

Hamiltonian (wie bei Pp4 3c ) und erhalt wieder ein recht komplexes Ergebnis, das sehr verkiirzt 
folgendermaBen lautet: 

2048 a 2 fr 63 (0 (l-2t) j-l+t) t e (l + 2t) 

PLS > 3 = 135(1 + 0" (4 ' 182) 

| 2048 a 2 F 366 (t) (-1 + t 6 (-1 + 2t) (l + 2t) 
135 (1 + t) 11 

1024 a 2 $ 23 (t) (-1 + t) t 6 (-1 + 2 t) (1 + 2 1) 
135 (1 + i) 11 

+f( 1 ,F 2 (t),F(t)M(t)). 

Daraus ist ersichtlich. dafi auch Pl s.3 wieder einen Beitrag zum numerischen Antcil beisteuert. Fiir -Fl-s^ 
erhalten wir 

„ o a 2 ln(a) a 2 ln(e) 

Fl-s,3 = -0.0018158 a 2 + -±-L -A2 + N L . S , 3 . (4.183) 

54 lUo 

Als Ergebnis fiir F^.s (Summc von Fl.s,i, F^.g^ und -Fl-s,3) erhalten wir 

F L . S = 0.0189735 a 2 + _ + ^ (41g4) 

18 36 



61 



4.2.8 Der numerisch auszuwertende Anteil 



Der verbleibende numerische Anteil aus Ns$,s, Nl.s,3 und N p i^ c ergibt sich zu 

N = Ns^ 5 + N L . S ,3 + N p ^ 3c (4.185) 



1 ( 512a 2 * 21 (t)t 256a 2 F 122 (t)t 



27 (1 - t) (1 + t) 5 27 (-1 + t) (1 + tf 



256 a 2 $21 1024a 2 F 366 (i) (l-2t) (-1 + t) 2 t (l + 2i) 



27 (-1 + *) (1 + i) 5 27 (1 + i) 10 

1024a 2 $ 63 ft) (1 — 2t) (~l + t) 2 t (1 + 2t) 
+ 27(l + <) 10 

512a 2 $ 23 (£) (-l + t) 2 t (-1 + 2t) (l + 2t) 
+ 27(l + t) 10 

Die numerische Integration wird mit der Gausschen Quadratur ausgefiihrt. Dabei erweist es sich als sehr 
hilfrcich, dafi der Integrand fur t — > und t — > 1 verschwindet. Dies wurdc durch das Hcrausnchmcn der 
k = 0- und fc = 1-Tcrme in den speziellen Funktionen erreicht. Das Endergebnis fur N lautet: 

N = 0.0030113 a 2 . (4.186) 



4.2.9 Ergebnis des Niedrigenergie-Anteils 



In alien Teilergebnissen zum Niedrigenergie- Anteil hatten wir der Einfachhcit halbcr Za durch a ersetzt. 
Wenn wir nun alle Teilergebisse aufaddieren und im Endergebnis a durch Za ersetzen, so ergibt sich: 



Ft. 



F, 



F n 



F„ 



F r xa + Fk-r 



(4.187) 



-F SE + F H + F s + F P 4 + F L . S + N 



103 

0.0400223(1) - 0.79565(1) (Za) 2 + — (Za) 2 ln((Za)~ 2 ) 

180 



2 (Za 2 103 , 2 . , , 

H — — H (Za) 2 ln(e). 

9e 180 V ' y ' 

Der Agi-Koeffizient des ln((Za) _2 )-Terms stimmt gerade mit dem bekannten Ergebnis von 103/180 
uberein [6]. Wir erinnern uns an das Ergebnis (4.73) fur den Hochenergie- Anteil 



F H (2P-_ 



1/2) 



6 



(Za) 



4177 103, 103, . . 2 

In 2 In (e) 

21600 180 w 180 w 9e 



und geben noch cinmal die Formcl fiir den Niedrigenergie-Anteil an. 

103 

F L (2P 1/2 ) = 0.0400223(1) - 0.79565(1) (Za) 2 + — (Za) 2 ln((Za)" 2 ) 

180 



2(Za) 2 103 , 2 , / v 
+A 97 L + 180 {Za) ~ ln(£) - 

In der Summc heben sich die Divergenzen in e hcraus, und wir crhaltcn als Summe der Hoch- und 
Niedrigenergie- Anteilc 



F(2P 1/2 



-0.126644(1) + (ZaY 



103 

-0.99891(1)+ — hi ((Za) 
180 



(4.188) 
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Die Kocfnzicntcn sind somit bestimmt zu 



A 40 (2P 1/2 ) 



G 



1 



\\nk {2P) = -0.126644(1) 



(4.189) 



Am 



0.99891(1), 



103 
180 



Die Werte fur und Aqi stimmen mit den bekannten Rcsultatcn iiberein (siehe [6]). Der Wert des 
Bethe-Logarithmus lnfc fur 2P-Zustande ist lnfc (2P) = -0.0300167089(3). 

Die numerischcn Integrationen crfolgen jewcils mit eincr Unsicherheit von 1 in der letzten Dezimalstcllc 
(d. h. in der 6. Nachkommastelle) . Im ganzen werden 10 numerische Integrationen ausgefuhrt. Wir geben 
daher die Genauigkeit des Ago-Koeffizienten mit Aqq = —0.99891(1) an (d. h. mit. einer Unsicherheit von 
1 in der 5. Nachkommastelle). 

Anmerkung: Es wird vermutet, daB die Integration sehr viel besser konvergiert als hier angegeben. 
Zum Bcispiel wurden bei der Integration des nichtrclativistischcn Dipol-Antcils allc 9 bekannten 
Nachkommastcllen des Bethe-Logarithmus in Ubcrcinstimmung mit dcm bekannten Resultat crhaltcn. Da 
aber beim ^460-KoefBzicnten die 5. Nachkommastelle bereits einer Genauigkeit von unter 1 Hz cntspricht, 
geben wir uns mit einer sehr vorsichtigen Abschatzung zufrieden. 
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Kapitel 5 

Die Selbstenergie des 2P-Zustands 
(j=3/2) zur sechsten Ordnung in 
Z_alpha 



5.1 Der Hoch- und Niedrigenergie-Anteil 

Die Rcchnungcn des vorangegangenen Kapitels konnen fiir das 2P 3 / 2 -Niveau in ahnlicher Weise wiederholt 
werden. Fiir den Hochencrgic-Antcil wird erhaltcn 



Fh(2P 3/2 



1 

12 



(Za) 



6577 29 



21600 90 



ln(2) 



29 
90 



ln(e) 



2 

97 



(5.1) 



Weil bcim 2P 3 / 2 -Niveau der Spin des Elektrons mit seincm Bahndrchimpuls zum Gcsamtdrchimpuls 
j = 3/2 koppelt, andcrn sich die Beitrage Fl-s, F rxa und Fk. r . Weil die relativistische Korrektur zur 
Energie fiir 2P 3 / 2 anders ist als fiir 2P 1 / 2 , andert sich F$e- Da der relativistische Hamiltonian SH die 
a-L-Kopplung enthalt, andert sich auch die relativistische Korrektur zur Wellcnfunktion S(f>. Damit erhalt 
auch Fsrj, einen anderen Wert. Die anderen Beitrage andern sich nicht. 

Es soil hier nicht auf die Rcchnungcn cingegangen werden. Im wescntlichcn wird dem bereits besprochenen 
2J\/ 2 -Zustand gefolgt. Der Niedrigenergie-Anteil ergibt sich zu 



Fl(2P 3/2 ) = --\nk a + {Zaf 



-0.58451(1) 



29 
90 



ln((Za)- 



29 
90 



ln(e) 



2 

97 



(5.2) 



5.2 Endergebnis fiir den 2P-Zustand (j = 3/2) 



In der Addition der Hoch- und Nicdrigencrgic-Antcile heben sich die 1/e- und ln(e)-Divcrgcnzcn gcradc 
hcraus, und wir erhalten als Ergcbnis fiir F 



F(2P 3/2 ) = 0.123356(1) + (Za)' 



29 

-0.50337(1) + — ln((Za)" 2 ) 
90 



(5.3) 



mit A40 = 0.123356(1) = 1/12 — 4/31nfco(2P) (hifco ist der Bethc-Logarithmus). Durch die Rechnung ist 
der Ago-Kocffizicnt bestimmt zu 

A 60 (2P 3 /2) =-0.50337(1). (5.4) 
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Kapitel 6 

Ergebnisse und Diskussion 



6.1 Zusammenfassung der Ergebnisse 



Es wurde in dicscr Arbeit die efw-Methode entwickelt, welche zur Bercchnung der (Ein-Schlcifcn-) 
Selbstenergie eines gebundenen Elektrons in hoherer Ordnung dient. Diese Methode wurde auf die 2Pi/ 2 - 
und 2P 3 / 2 -Zustande angewandt. 

Fiir die F-Faktoren der Selbstenergie (zur Definition von F siehe Kapitel 1.3) erhalten wir 



und 



F(2P 1/2 ) = ~ - | lnfc (2P) + (Zaf 



P(2P 3/2 ) - -L -ibxko(2P) + (Zaf 



1 03 

-0.99891(1)+— In ((Za)~ 2 ) 



29 

-0.50337(1) + — In ((Zq)" 2 ) 



(6.1) 



(6.2) 



wobei der Bethe-Logarithmus gegeben ist durch lnfc (2P) = —0.0300167089(3). Die A± - und Agi- 
Kocffizicntcn stimmen mit den bckannten Rcsultaten iiberein (diese Kocffizicntcn sind in [C] fiir P und 
S Zustande aufgefiihrt, fiir die P Zustande siehe Kapitel 1.3). Die ^go-Kocffizicntcn sind bestimmt zu 



^6o(2Pi/ 2 ) = -0.99891(1) 



(6.3) 



und 



A> (2P 3/2 ) = -0.50337(1). (6.4) 

Diese Resultate sind konsistent mit den aus der Extrapolation von numerischen Rcsultaten gewonnenen 
Werte fiir die Ago-Koeffizienten [14], 

A 6 o(2Pi /2 ) = -0.8(3) (6.5) 

und 

A 6 o (2P 3/2 ) = -0.5(3). (6.6) 

Die Bestimmung der Ago-Koeffizienten mit verbesserter Genauigkeit war neben der Entwicklung der efw- 
Methode das Hauptziel dieser Arbeit. Mit der gcnaucn Kcnntnis der Kocffizicntcn wird cs moglich, neue 
(genauere) Werte fiir die Lamb-Shift der beiden 2P Zustande anzugeben (Kapitel 6.4). 



6.2 Uberpriifung der Rechnung 

Lamb-Shift-Rechnungen allgemein sind sehr umfangreich und die Ausdriicke in den Zwischenschritten 
recht grofi. In dieser Arbeit iiberstiegen die Integranden im Hochenergie-Anteil bei den Matrix-Elementc 
eine Lange von 1000 Termen, im Niedrigenergie-Anteil traten Ausdriicke vergleichbarer Lange auf, 
die dann vermittels der Zusammenhangsformcln vcrcinfacht wurden. Ein einziger Rechenfehler kann 
das Ergebnis verfalschcn. Dahcr stellt sich die Frage nach den Uberprufungsmoglichkeiten. Bei der 
vorliegcnden Arbeit wurden folgcnde Prufungen vorgenommen: 
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1. Die Matrix-Elemente fur den Hochenergie-Anteil wurden auf zwei verschiedene Arten 
(computergestiitzt und von Hand) ausgewertet. 

2. Im Hochenergie-Anteil sind der Kocffizicnt der ln(e)-Divcrgenz und der 1/e-Divcrgcnz bckannt und 
konnen iibcrpriift werden. 

3. Im Niedrigenergie-Anteil kann man fiir jeden einzelnen Beitrag die ln(e)-Divergenz berechnen und 
iiberpriifen. Dies wurde fiir jeden einzelnen Beitrag getan. Die Summe der Divergenzen mufi sich 
gerade mit dem Hochenergie-Anteil herausheben. 

4. Man kann zeigen, dafi das Matrix-Element 
mit 

P*>) = Wa^?- 1 -^e- iS -» 

fiir ui — > verschwindet. Dahcr mufi auch die Summe der einzelnen Beitrage des Nicdrigcncrgic- 
Antcils fiir w — > (t — ► 1) vcrschwinden. Nach der FW- Transformation treten jedoch noch 
vcrstccktc Beitrage auf, da die o'-Matrizcn obere und unteren Komponente der Spinoren 
miteinander verbinden. Wenn man die verstccktcn Beitrage jedoch mitcinbezicht, so mufi sich 
der konsistcnte Wert Null fiir die Summe ergeben. Man kann den Wert der Beitrage fiir t — ► 1 
sogar einzeln berechnen, wodurch eine zusatzliche Uberpriifung des Resultats moglich ist. Diese 
Priifungsmoglichkeit wurd ebenfalls in die Tat umgesetzt. 

5. Alle im Nicdrigcncrgic-Anteil verwandten Summations-Rcgcln wurden an gceignetcn Beispiclcn 
gepriift. 

6. Der nichtrelativistische Dipol-Anteil wurde auf zwei verschiedene Arten (im Impuls- und Ortsraum) 
berechnct. 



Diese Priifungen beziehen sich auf die Korrckthcit der Rechcnschrittc. Es ergibt sich noch ein weiteres 
Problem, denn auch die verwandten Computcr-Algcbra-Systcme sind (wahrscheinlich) nicht frci von 
Inkonsistenzcn. Die Uberpriifungsmoglichkcitcn in Bezug auf diese Systcme sind eingeschrankt. Man 
kann allcnfalls die cigencn Programme so gcstalten, dafi die einzelnen Rechcnschrittc schr gut kontrolliert 
werden. 



6.3 Abschatzung hoherer Koeffizienten 

P. Mohr hat F fiir grofie Z aus numerischen Rechnungen gewonnen [14]. Fiir Z = 5 ist die Konvergenz 
bereits sehr langsam, fiir Z < 5 werden dort keine Resultate angegeben. Aufgrund der nichtanalytischen 
Terme in a ist eine Extrapolation der Ergebnisse in den Bereich fiir kleine Z nur eingeschrankt 
moglich. Kennt man aber die Koeffizienten von F bis zu einer bestimmten Ordnung in a, so lassen 
sich durch Angleichen einer Funktion, die die nachsthoheren Koeffizienten enthalt, Abschatzungen fiir 
diese Koeffizienten gewinnen. Wir gleichen dazu die Funktion 

nl7 , HZ) - A 40 - {Zaf (A m + Mi In [(Za)- 2 ] ) 

G{Z} (Zaf _ (6J) 

A n + (Za) [a 80 + A 81 In [(Za)- 2 ] + A S2 (in [(Za)" 2 ]) 2 ] + 0((Zaf) 

an die numerischen Ergebnisse an (mit den in dieser Arbeit bestimmten, neuen Werten von Aqq). Daraus 
ergeben sich die folgenden Abschatzungen fiir die nachsthoheren Koeffizienten. 

49 

A 70 (2P l/2 ) a -7T, (6.8) 
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Die Angabe dieses Koeffizienten erfolgt in dieser Form, da gezeigt werden kann, dafi A ungera( jo.o gleich 
einer rationalen Zahl mal 7r) ist. 

A 8 o(2Pi/ 2 ) « 0.1, (6.9) 



Fur den 2P 3 / 2 -Zustand ist 



A 81 (2P 1/2 ) « -1.5, (6.10) 

A 82 (2P 1/2 ) w 0.2. (6.11) 
39 

A 7 o(2P 3/2 ) « (6-12) 

^8o(2P 3/2 ) « -1.7, (6.13) 

A 8 i(2P 3/2 ) « -0.6, (6.14) 

A 82 (2P 3/2 ) » -0.2. (6.15) 

Den Fehler dieser Abschatzungen sollte man aufgrund der systematischen Unsicherheit des Verfahrens 
mit mindestens der Halfte des Absolutwertes in jedem einzelnen Fall ansetzen. 

Fiir die Werte von G\ = G(Z = 1) ergeben sich die folgenden Abschatzungen: 

Gx(2P 1/2 ) = 3±2, (6.16) 

und 

Gi(2P 3/2 )=2±l. (6.17) 

Diese Werte von G werden fiir die Angabe des theoretischen Fehlers des Selbstenergiebeitrages bei der 
Berechnung der Lamb-Verschiebung benotigt. 

6.4 Neue Werte fur die Lamb-Verschiebung 

Fiir die Lamb-Shift 5pLamb verwenden wir die Definition 

2 

E = M + m r [f(n,j) - 1} - ^ [f(n,j) - if + 6E hamb + E Ms . (6.18) 

E ist die Encrgie des Zwcikorpersystems und f(n,j) die dimensionslosc Dirac Energic. Fiir die in 
dieser Arbeit behandcltc Ein-Schleifcn-Sclbstenergic muB zur Bcstimmung der Lamb-Verschiebung die 
Abhangigkeit der einzelnen Koeffizienten von der reduzicrtcn Masse des Systems miteinbezogen werden. 
Diese Abhangigkeit ist bekannt [6]. Aufier der Ein-Schlcifcn-Sclbstenergie tragen folgcndc Effcktc zur 
Lamb-Verschiebung bei : 



1. Zwei-Schleifen-Korrekturen zur Sclbstcnergic. Der Hauptbeitrag ist von der Ordnung a(Za) 4 . Es 
ist auch ein logarithmischcr Bcitrag der Ordnung a 2 (Za) 6 ln((Za)~ 2 ) bekannt [5]. 

2. Drei-Schleifen-Beitrag zur Selbstenergie. Der Hauptbeitrag dieser Korrektur ergibt sich durch die 
Drei-Schlcifen-Beitrage zum anomalen magnetischen Moment des Elektrons [5]. Der Effekt skaliert 
mit a 3 (Za) 4 . 

3. Vakuumpolarisation. Die Vakuumpolarisation geht fiir P-Zustande nur in hohcrer Ordnung 
(a{Za) 6 ) ein [6]. 

4. Eine mit (m r /M) 2 (Za) A skalierende Ruckstofikorrcktur, die aus dem Brcit-Hamiltonian folgt und 
eine leichte /-Abhangigkeit der Energieniveaus bedingt [6]. 

5. Die Salpeter Korrektur (relativistischc Ruckstofikorrcktur) in {Za) 5 [6]. 
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6. Relativistische Rvickstofikorrekturen hoherer Ordnung ((Zaf m r /M) [24]. 

Durch Abschatzen hdhcrcr Tcrmc in a kann fur die cinzclncn Beitrage ein "theoretischer Fehler' 
angegeben werden. Fiir den 2P 1 / 2 -Zustand ergibt sich folgende Tabelle, 



Lamb-Vcrschiebung des 2P 1 / 2 -Zustands 


Beit rag 


\ \ /Avf i n 1 r I— 1 n 

wen m kxiz 


Fehler in kHz 


T^ji Ti- Scb 1 oi foii-Sol bst oil ori^i o 

1 Jill k w 7\_.±J.lv^ll.l_'i.l X L/O L cLIVj 


-12846 Q 


0.1 


Zwei-Schleifen-Selbstencrgic 


26.0 


0.1 


Drci-Schlcifcn-Sclbstcncrgic 


-0.2 


vcrnachlassigbar 


Vakuumpolarisation 


-0.3 


vcrnachlassigbar 


(Za) 4 Riickstofi 


2.2 


vcrnachlassigbar 


(Zaf Riickstofi (Salpeter) 


-17.1 


0.1 


{Zaf Riickstofi 


0.4 


vcrnachlassigbar 


Summc fiir 2P 1 / 2 


-12836.0 


0.3 



und fiir 2P 3 / 2 haben wir 



Lamb-Vcrschiebung des 2P 3 / 2 -Zustands 


Bcitrag 


Wert in kHz 


Fehler in kHz 


Ein- Schlcifcn- Sclbst cnergie 


12548.0 


0.1 


Zwei-Schleifen-Selbstencrgic 


-12.8 


0.1 


Drci-Schlcifcn-Sclbstcnergic 


0.1 


vcrnachlassigbar 


Vakuumpolarisation 


-0.1 


vcrnachlassigbar 


(Za) 4 Riickstofi 


-1.1 


vcrnachlassigbar 


(Zaf Riickstofi (Salpeter) 


-17.1 


0.1 


(Zaf Riickstofi 


0.4 


vcrnachlassigbar 


Summc fiir 2_P 3 / 2 


12517.5 


0.3 



Bei den theoretischen Unsichcrheiten addieren wir die Absolutwerte der Fehler, da sich theoretischc 
Unsichcrhcitcn nicht wie statistisch verteiltc Mefiwertc verhaltcn. Die ncuen Wcrtc fiir die Lamb- 
Verschiebungen der 2P-Zustande sind 

(5P Lamb (2P 1/2 ) = -12836.0(3) kHz (6.19) 

sowie 

<5P Lamb (2P 3/2 ) - 12517.5(3) kHz. (6.20) 

Der theoretische Fehler fiir beide Endergebnisse liegt bei unter 300 Hz. Der neue theoretische Wert fiir 
die Feinstruktur ergibt sich aus Gl. 6.18 zu 

SE is = P(2P 3/2 ) - P(2P 1/2 ) = 10969043.2(6) kHz, (6.21) 

wiedcrum mit der Mafigabe, die Absolutwerte der einzelnen theoretischen Fehler zu addieren. 

6.5 Diskussion der epsilon-fw-Methode 

Die in dicscr Arbeit vorgestclltc efw-Mcthodc dicnt zur Bercchnung der Selbstcncrgic eincs 
gebundenen Zustands in hoherer Ordnung. Die Rechnung wird in einen Hoch- und einen 
Niedrigenergie- Anted aufgespalten. Der Hochenergie- Anted wird durch Einsatz von Computer-Algebra 
vereinheitlicht behandelt. Im Niedrigenergie- Anted hat die efw-Methode aufgrund der Foldy-Wouthuysen- 
Transformation folgende Vorteile: 

1. Klare Trcnnung des Hauptbeitrages von den relativistischen Korrekturen, 

2. Einsatz des nichtrelativistischen Propagators, 
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3. Aufteilung der Rechnung in Teilbeitrage. 

Ultimativcs Ziel der computcrgcstiitztcn Rcchnungcn ist automatisierte Berechnung von 
quantcnclcktrodynamischen Korrckturcn. Nur durch cine solchc Automatisicrung der Rechnung 
ist es moglich, Korrekturen hoherer Ordnung zuverlassig zu berechnen (wegen der Lange der Terme). 
Von besondcrcr Wichtigkcit blcibt jedoch die Wahl einer adaquaten Rcchenmcthode. Bei Wahl einer 
falschen Methode "explodieren" die Ausdriicke in mehrere zehntausend Terme, die auch mit schnellen 
Rechnern nicht mehr handhabbar sind. Daher ist es unter anderem wiinschenswert, die Rechnung in 
Teilschritte aufzuspalten. 

Wie wir am Beispiel der Integrationsroutinen zur Berechnung der Matrix-Elemente (siehe Abschnitt 
4.1.1) sehen, kann durch geschicktes Programmieren haufig die Rechenzeit um mehrere Groficnordnungen 
vcrringcrt werden. Die entwickcltcn Integrationsroutinen zur Berechnung von Matrix-Elcmcnten kommen 
mit sehr wenigen Rcgeln aus und sind daher schncll. Die in Mathematica eingebautcn Routincn decken 
zwar eine groficre Klassc von Integralen ab, sind aber fur den spezicllcn Einsatz viel zu langsam. 



09 



Kapitel 7 
Ausblick 



7.1 Anmerkung zum Vergleich mit dem Experiment 

Zum direkten Vergleich von Theorie und Experiment sollte auf dem kHz-Niveau noch folgender Effekt 
miteinbezogen werden, fur den hier dcr Name hfs-fs-Kopphmg vorgeschlagen wird. Der Effekt wird 
beispielsweise in [23] fur myonischen Wasserstoff behandelt. Er beeinfluBt die relative Lage von Zustanden 
mit denselben n, I, und F Quantcnzahlcn, die aber verschiedene j-Quantcnzahlen aufweisen. Der Effekt 
wird bereits in dem klassischen Buch von Bethe und Salpeter [9] beschrieben. Durch den Effekt wird der 
2-p3/2j F = 1 Zustand energetisch abgesenkt, der 2P 1 / 2 , F — 1 Zustand jedoch cncrgctisch angehoben 
um jeweils 2.5 kHz. Dcr Effekt beeinfluBt nicht die relative Lage der 2P 3 / 2 , F — 2 und 2P 1 / 2 , F = 
Zustande. 

7.2 Hochprazisions-Tests der QED 

Durch den neuen theoretischen Wert fur den 2P 1 / 2 Zustand wird cs moglich, von der 2P Lamb- 
Verschicbung, die als Refcrenzpunkt dicnt, die 2S Lamb-Vcrschiebung sehr genau zu bestimmcn (Mcssung 
von Hagley und Pipkin, [13]). Andererseits wird die 2«S , -Lamb-Vcrschicbung im Frcquenzkcttcnlabor 
durch Vergleich der nS'-Ubergangsfrequenzen ermittelt. Dadurch ist ein direkter Vergleich von zwei vollig 
unterschiedlichen experimentellen Methoden moglich. 

Allgemein gilt fur die hier vorgestellte efw-Methode, dafi analytische Rechnungcn zur Selbstenergie 
in hSherer Ordnung sehr vereinfacht werden.. In leicht abgewandelter Form sollte die Methode zur 
Berechnung von Strahlungskorrekturen zur Hyperfeinstrukter in Myonium und Positronium sowie zur 
Selbstenergie in Helium einsetzbar sein. 

Weitere Anwendungsmoglichkcitcn crgcben sich fur P-Zustandc mit hohcrcn Hauptquantcnzahlcn (n = 3 
und n = 4). Die Computerunterstiitzung ist hier sehr hilfrcich, da die aufgrund hohcrcr Hauptquantcnzahl 
langcren Ausdriickc fur die Wcllcnfunktioncn problcmlos behandelt werden konncn. 

Nach Wissen des Autors existieren fur Zustande mit n > 3 zur Zeit weder numerische noch symbolischc 
Rechnungcn in hohcrcr Ordnung. Durch Rechnungen in hoherer Ordnung konnte die Lamb-Shift dieser 
Zustande und damit die Feinstruktur sehr genau bestimmt werden. 
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